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ABSTRACT 

We consider a particular class of threshold corrections to Yukawa couplings and mass rela- 
tions in the MSSM and supersymmetric grand unified models. We give a complete treatment 
of Yukawa coupling thresholds at the unification scale M GUT and the effective supersymmetry 
scale M SUSY and apply them to corrections to the tree-level prediction yb(M GUT ) — y T {M GUT ) in 
minimal supersymmetric SU(5). We apply both gauge and Yukawa coupling thresholds to gauge 
unification and the above Yukawa unification condition to find predictions for the top quark mass, 
Mt, the superheavy vector boson mass M v and the colored Higgs triplet mass M Hs . We discuss 
the dependencies of M GUT and M H3 on a s (Mz), M t and the sparticle spectrum as well as those 
of M t on tan/3, a s (Mz) and the bottom quark mass, Mb- The effect of the Yukawa coupling 
thresholds on M t are given for representative sparticle spectra. We describe the quantitative 
differences between these effects for low and high tan/3. We also give new bounds on superheavy 
masses, incorporating proton decay as well as unification constraints, the former leading to a 
lower bound on a s . 
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1. Introduction 



In light of the current renewal of interest in renormalization group (RG) con- 
straints arising from supersymmetric grand unified theories (SUSY-GUTs), it has 
become increasingly important to quantify corrections to GUT scale predictions. 
Recently several two loop renormalization group analyses of SUSY-GUTs with 
soft supersymmetry-breaking induced via minimal N=l supergravity have been 
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performed. In certain scenarios, these analyses have made predictions for the 
sparticle spectrum and low energy parameters arising from GUT scale constraints. 
Emphasis should be placed on determining the uncertainties in these predictions 
and their dependence on the unknown heavy mass spectrum as well as the details 
of supersymmetry breaking. However, only some of these analyses treat the one 
loop threshold matching conditions which are necessary for a consistent two loop 
analysis of the gauge coupling unification and few consider thresholds occuring in 
Yukawa couplings and mass parameters in general. Usually threshold effects are 
considered only for the gauge couplings. 1 We find that the Yukawa coupling thresh- 
olds are typically at least as important as the gauge coupling threshold effects. 

The purpose of this article is to present a complete treatment of Yukawa thresh- 
old corrections in the Standard Model (SM), the Minimal Supersymmetric Stan- 
dard Model (MSSM) and the minimal supersymmetric SU(5) GUT. We also review 
the gauge coupling threshold corrections in these models and attempt to quantify 
both types of correction terms. These effects are applied in the context of gauge 
and third generation Yukawa unification in minimal SU(5), and are used to probe 
the superheavy spectrum and the top quark mass. For definiteness, we consider 
the scenario in which these unifications, modulo grand unified scale thresholds, are 
used to solve for the physical top mass M t , the colored Higgs triplet mass M Hz 
and an effective GUT scale M GUT . These solutions are considered for experimental 
ranges of values of sin#vf/> the strong coupling a s and the physical bottom quark 

1 For a discussion of this indirect influence of gauge coupling thresholds on the prediction of 
m b{Mz) from Yukawa unification, see Ref. 8. 
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mass M;, and central values of the remaining weak scale parameters. We also give, 
especially for the M t solution, the dependence on tan/5, the ratio of the VEVs 
of the two Higgs doublets in the MSSM (see Appendix A for conventions). We 
highlight the importance of threshold effects in these solutions by displaying their 
quantitative dependence on the sparticle and GUT scale spectrum. 

We find significant variation of the solutions for M Hs and M GUT with top masses 
in the range 130 < M t < 200 GeV. The corresponding weak threshold corrections 
are between 1 and 40% for the M GUT solution and between 6 and —95% for the M Hz 
solution. The latter solution is more sensitive to the large a s error bars, two loop 
effects and boundary corrections at Mz and M SUSY , particularly since M H3 itself 
only appears in threshold matching functions. The decline of the M Hs solution 
with Mt is to be contrasted with the increase in the lower bound on M Hs from 
proton decay in which a large top Yukawa enhances the corresponding effective 
dimension 5 operator through renormalization group effects. We allude to the 
possibilities for constraining the minimal SUSY-SU(5) model using proton decay 
limits without introducing theoretical prejudices about the GUT scale spectrum. 
We also re-emphasize the sensitivities of M GUT and M Hz on the gluino and Higgsino 
masses respectively. 

The effect of the Yukawa thresholds is seen in the Mt solution or, equivalently, 
m b(Mz) f° r fixed Mt- We discuss the robustness of the former solution compared 
to the latter in the context of both Yukawa threshold corrections at M SUSY and 
M GUT . We find that these GUT scale corrections to Mt typically range from +10 
to —5% for extreme ranges of the ratio of the superheavy vector and adjoint scalar 
masses. We dissect the various supersymmetric Yukawa coupling threshold correc- 
tions, in particular the various sources of enhancements in the large tan (5 region 
coming from certain vertex diagrams involving gluinos and charginos. Quantitative 
effects for M t are given for sample sparticle spectra for low, intermediate and large 
tan p. Characteristically, we find positive corrections to M t of order 2-7% for low 
tan/3 = 1.5 and for both signs of the supersymmetric Higgs mass parameter \x H . 
For higher tan/3 and ji H > 0, we find typically negative corrections of up to 10% 
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for tan/3 = 15. For tan/3 <; 40, the corrections are larger and highly spectrum 
dependent. The inclusion of the Yukawa coupling thresholds can also significantly 
modify perturbativity constraints on third generation Yukawa couplings. 

The importance of Yukawa coupling thresholds has come to the attention of 
several authors recently.' 9 111 Loop corrections to low energy mass parameters have 
been treated in certain cases such as the well-known case of radiative corrections 

[12] 

to the Higgs boson masses in the MSSM although usually not explicitly in the 
form of threshold conditions on running mass parameters except in the case of 
the Standard Model' 131 and for y^ and y T in ordinary GUTs.' 141 However, in the 
context of renormalization group analyses of high scale predictions they have not 
been treated generally. 

This work is organized as follows. We first consider in Section 2 the gen- 
eral treatment of Yukawa and mass parameter thresholds. This analysis provides 
a framework which can be applied to more complicated mass and mixing angle 
textures. In Section 3 we give a complete calculation of the GUT scale match- 
ing functions for the Yukawa couplings in minimal SUSY-SU(5). In the spirit of 
quantifying uncertainties in SUSY-GUT predictions, we consider in particular the 
one loop threshold corrections to the tree-level relation yb(M GUT ) = y T (M GUT ) and 
highlight its sensitivity to the top Yukawa coupling yt{M GUT ) and the possible split- 
ting of the colored Higgs triplet superfields from the other superheavy fields. This 
is one of many threshold corrections to various mass relations discussed recently 
in Refs. 15 - 19 when a specific GUT or string model is taken to generate certain 
textures in the fermion mass matrices at the unification scale. We also review the 
GUT scale matching functions for the gauge couplngs. 

In Section 4 we give the most important contributions to the Yukawa coupling 
matching conditions at M SUSY in the MSSM in the case of one light Higgs doublet 
in the effective SM below M SUSY . We show the origin of the large enhancements 
which occur for large tan (3. We discuss the sensitivity of these corrections to the 
soft supersymmetry breaking parameters of the model. A complete calculation of 
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these thresholds for a general sparticle spectrum is given in Appendix A. We also 
review the gauge coupling thresholds at M SUSY as well as the top mass dependent 
gauge and Yukawa thresholds when the top quark together with the light Higgs 
and W and Z bosons are integrated out at Mz- Finally we describe our procedure 
for extracting the heavy quark pole masses. 

In Section 5 we apply these threshold corrections to analytic predictions for 
M GUT = (M 2 M s )i, M H . A and m t (Mz), where M v and M E are the superheavy 
vector and the superheavy adjoint scalar masses, respectively. The first two pre- 
dictions are essentially determined by gauge coupling unification. Generally gauge 
unification can be used to determine the unification scale M GUT , the value of the 
couplings at M GUT , a G , and either a s (Mz), sin 9w(Mz) or M Hy In the latter case 
the mismatch of the gauge couplings at the GUT scale determines the colored Higgs 
triplet mass through threshold corrections. We choose this somewhat unconven- 
tional case as it is more useful in discussing proton decay constraints. We quantify 
the importance of threshold corrections at Mz and M SUSY on the M GUT and M Hz 
solutions. In particular, we describe the importance of the top mass dependence 
of sin 2 6w(Mz) and the gluino and Higgsino masses, respectively, on these two 
solutions. For the last prediction, we give a semi-analytic solution for mt which 
allows one to identify the relative importance of the various correction terms. We 
note the robustness of the mt solution in all but the extreme large tan/3 region 
and relate this to the attraction of the infrared quasi- fixed point in yt- For all the 
solutions we show their dependence on the experimental uncertainties in low scale 
parameters. 

In Section 6 we briefly review the proton decay bound on M H3 and indicate 
how conservative bounds can be improved when applied to the context considered 
here. We also review how perturbativity arguments are used to put theoretical 
bounds on the GUT scale spectrum. We apply these, together with the condition 
that all masses lie below the Planck mass M Planck , only to delineate the cases in 
which a perturbative analysis is appropriate. 
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We give the full two loop numerical solutions for M GUT , M„ 3 and M t in Section 
7. We give the extreme ranges of M GUT and M H . A consistent with gauge coupling 
unification and low energy data. We display the effects of the weak and supersym- 
metry scale threshold corrections as well as the error bars on a s of these solutions 
in graphical form. The effects of the thresholds at M SUSY and M GUT on the Mt 
solution are separately discussed. In the former case we quantify the effects in the 
low, intermediate and high tan (3 regions for sample sparticle spectra. In the latter 
case we show how the GUT scale spectrum can modify the prediction for Mt as 
a function of tan j3. We also show the dominant effects of the uncertainties in M5 
and a s on Mt. 

Finally we make two general remarks. The first concerns our restriction to 
minimal SUSY-SU(5). This restriction is chosen merely for its relative simplicity. 
The minimal case may also be indicative of what one may expect in more general 
GUTs. Thus we do not address the less appealing features of minimal SU(5) such 
as the doublet-triplet splitting problem and the unsuccessful mass relations for 
the first and second generation fermions. These difficulties can be resolved by 
complicating the GUT scale superpotential.' 201 Instead we regard minimal SUSY- 
SU(5) as a toy model for investigation. In general, one expects the threshold 
corrections in more complicated models to be at least as large as those in the 
minimal case. 

Second, we would like to advocate a particular philosophy for the treatment of 
threshold corrections in minimal subtraction (MS) schemes. From the Appelquist- 

[21] 

Carrazzone decoupling theorem we expect the physics at energies below a given 
mass scale to be independent of the particles with masses higher than this thresh- 
old. However, MS schemes are not physical in the sense that they are scale depen- 
dent and mass independent so that the decoupling theorem is not manifest. As 
described in Refs. 22 - 24, one implements the decoupling in MS schemes by for- 
mulating a low energy effective theory obtained by integrating out the heavy fields 
to one loop. The effect of this procedure is to give relations between renormal- 
ized parameters just above and below the particle mass(es), the so-called matching 
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functions, and to modify the various (3 function coefficients so that in the lower 
scale theory the contribution of the particle(s) to these coefficients is removed. 
The running parameters are typically discontinuous at the boundary at which the 
matching function is applied unless one tunes the boundary scale. In either case 
reliable values for the parameters of the theory are obtained asymptotically away 
from the boundary. 

In the treatment of gauge coupling thresholds, particularly in the MSSM and in 
GUTs, many authors have identified a threshold with each sparticle (often assum- 
ing small mixings or large degeneracies) or superheavy particle mass, and modified 
the beta function coefficients between each threshold. Here the couplings are usu- 
ally taken to be continuous at the boundary, which is justified since the gauge 
matching functions in supersymmetric theories in the dimensional reduction (DR) 
scheme vanish at the mass of the particle defining the boundary. However this pro- 
cedure can be complicated when the effect of integrating out the particle breaks the 
gauge symmetry of the higher scale theory. In the case of split multiplets (such as 
the Higgs 5 of SU(5)), one typically must use different sets of f3 functions depend- 
ing on different mass orderings of the spectrum. Also, when Yukawa thresholds 
are included the Yukawa and mass parameters are in general not continuous at 
threshold boundaries even in DR schemes. 

[23] 

Instead we use the arguably simpler scheme advocated by Hall in which 
one integrates out together all particles with similar masses at a single scale. If 
one uses one loop matching functions then this is justified as long as one loop (5 

[23] 

functions can be used between the different particle masses. This is generally 
the case. This is simpler both from the analytic and numerical standpoint when 
it is applied to the MSSM and its grand unified extensions. For example, in the 
analysis below in which one has a light SM Higgs doublet below M SUSY) the top 
quark as well as the W , Z, and SM Higgs are integrated out at Mz, all sparticles 
including the heavier Higgs doublet can be integrated out at a fixed scale, M SUSY 
(with complicated matching functions incorporating the details of the spectrum) 
and all superheavy particles at M GUT . With more generality and without a loss of 
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accuracy, the proliferation of scales and sets of (3 functions is reduced to at most 
three. 



2. Yukawa Thresholds: The General Case 

The treatment of Yukawa coupling and mass parameter thresholds in minimal 
subtraction (MS) schemes is analogous to the treatment of gauge thresholds given 

[221 [231 f24l 

by Weinberg, Hall and Ovrut and Schnitzer. For completeness we outline 
the procedure for Yukawa coupling and mass thresholds below. 

Consider a generic gauge theory in which the kinetic and Yukawa contributions 
to the bare Lagrangian are 

l~-kin = £gauge 

(2-1) 

where ipi is a generic fermion of type i, <p a is a generic scalar field with vacuum 
expectation value v a , and My = Y^v a - We assume that it is possible to decompose 
the fields ipi and (f> a into light components ipi, (f> a and heavy components and 
$ a When the heavy fields are integrated out to one loop we generate the following 
low energy effective Lagrangian, 

4£ = tfllge + *$Li + ^ Ri ^Rj + ^ a Zfd^ h + ... 

£y f = ~ ' 1,^1 i<.i ' ^LiV^Rjk + h-C. + NR , 

where NR indicates induced nonrenormalizable interactions. The matrix parame- 
ters Z^ L R ^ can be written as 

Z {LtR) = 1 + K {LiR) , (2.3) 
where K^ L R ^ come from light fermion wavefunction renormalization diagrams with 
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heavy fields in the loop. The parameters in 

= 1 + , (2.4) 

comes from analogous light scalar wavefunction renormalization. The matrices M 
and Y can be written as 

M = M + 5M = Mil + K M ) , 

(2.5) 

Y = Y + 5Y = Y(l + K Y ) , 

where SM and SY come from one loop light field self energy and three point Yukawa 
vertex diagrams, respectively, with heavy fields in the loop. 

To properly normalize the kinetic terms one must first diagonalize the hermi- 
tian matrices Z( LR ^ and then rescale each bare fermion field by the appropriate 
eigenvalue. Let = \j\ j ZfU L where Z^ is diagonal and Ul is unitary, with similar 
definitions for U R and U^. We redefine the fields via 

Via = (Z d Li )-U^ Lj , 

Vm-iZ^U^Rj, (2.6) 

In terms of the primed fields, the kinetic terms are canonical and the bare effective 
mass and Yukawa matrices are 

M eff = {Z d L )-?U L {M + 5M)U R {Z d R y? , (2.7) 

yeffa = {Z d L )^U L (Y + SY) b U R (Z d R )^ul ba (Z d a )^ , (2.8) 

To get the one loop matching functions one must turn these bare relations into 
renormalized ones. The divergent parts in M e ^ and Y e ^ a are just the difference 
between the self energy counterterms of the effective theory and those of the full 
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theory. So by replacing all parameters in (2.7) and (2.8) by their finite parts 
(indicated by a bar over the quantity) as defined in MS and neglecting higher 
order corrections, one obtains the renormalized relations: 

Mf(v) = {U L M^)U% + My(jjL)(K Mij - \{K d Li + Kj)) , (2-9) 

Yg fa (li) = (U L Y%)U^ul ba + Y^){Tt Yij - \{^ Ll + K d Rj + tfjj) , (2.10) 

In practice M e ^ must be diagonalized by a biunitary transformation even if the 
original theory is written in terms of mass eigenstates. In the case of the MSSM, 
for example, this leads to threshold corrections in the renormalized Cabibbo- 
Kobayashi-Maskawa matrix at the effective scale M SUSY when the superparticles 
are integrated out. The effect of such thresholds, of potential relevance to RG 

[25] 

analyses of GUT scale mass textures, will be discussed in a future work. 



3. Yukawa and Gauge Thresholds in Minimal SUSY-SU(5) 

Next we apply the preceding formalism to the one loop GUT threshold correc- 
tions to the tree-level relation yb{M GUT ) = y T {M GUT ) in minimal supersymmetric 
SU(5). This is of particular interest since this condition strongly constrains the 
allowed range of the top quark mass in an experimentally accessible region.' 26 291 
Here we find two contributions, one proportional to the square of the top Yukawa 
coupling coming from integrating out heavy color triplet Higgs and Higgsino fields, 
the other proportional to the GUT scale gauge coupling coming from integrat- 
ing out superheavy vector bosons. Since the soft supersymmetry breaking mass 
parameters are much smaller than typical GUT scale masses we can work in an ap- 
proximately supersymmetric formalism in which both light fields along with their 
superpartners are treated as massless in loops involving superheavy fields. We 
therefore use supergraph methods to simplify the calculations. 
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We start with a superpotential of the form 



P = ^ a Y^ X ab Hl l) - \e abcdeX ah Y^ X cd H^ 



+ M 2 H^ a H^ + X 2 H^ b X 1) + ^ Tr S 3 + ^ Tr S 2 , 



Mi 



(3.1) 



where iJ^, .ff( 2 ), x and £ are SU(5) superfields transforming as the 5, 5, 5, 10 
and 24 dimensional representations, respectively. In the usual way we associate 
these superfields with the SU(3) x SU(2) x U(l) decomposition by 



X 



1 



/ 







-U C 2 


-U 1 


-D'\ 









m 


-u 2 


—D 2 













—D 3 




u 1 


u 2 


f/ 3 





-E c 


V 


D 1 


D 2 


D 3 


E c 


o ) 



(3.2) 



and ^ T = (D^,E,— v), where lowercase Greek/Latin indices denote SU(3) and 
SU(5) indices, respectively. In component fields the charge conjugate Weyl fermion 
fields with be related to the right-handed component of Dirac fermions. The five 
dimensional Higgs supermultiplets are likewise split into color triplet and SU(2) 
doublet components according to 

« = (Hg>, H^-, and = (H^ a , H^+, H&) , (3.3) 



where the 3 denotes the strongly interacting triplets. These superfields couple to 
the SU(5) vector supermultiplet in the usual way. In addition to the component 
field interactions obtained from (3.1), one also can add the most general set of SU(5) 
invariant soft supersymmetry breaking terms. In the threshold calculations that 
follow, however supersymmetry as well as electroweak symmetry breaking effects 
are unimportant as the corresponding mass scales are much smaller than M GUT . 
The SU(5) symmetry is broken without breaking supersymmetry when the adjoint 
Higgs field gets the vacuum expectation value (£) = Vdiag(l, 1, 1, — 3/2, — 3/2), 
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where V = 2M\j\\. This gives the following mass spectrum in the minimal model: 
two degenerate Higgs triplet superfields of mass M H3 , degenerate X and Y lepto- 
quark gauge superfields of mass M V) color octet and SU(2) triplet superfields of 
mass M E along with a singlet of mass 0.2M E . The electroweak doublets remain 
light as long as the fine-tuning constraint X1M2 — 3A2M1 m M EW is satisfied, where 
M EW is a typical electroweak scale mass. This can of course be relaxed in the miss- 
ing doublet model for example,' 201 but we limit our discussion to the minimal model 
for simplicity. 

The GUT scale threshold corrections are obtained by integrating out the su- 
perheavy X and Y gauge supermultiplets, the Higgs triplet, and the Higgs ad- 
joint superfields to obtain an effective MSSM. These were obtained for the gauge 
couplings long ago in Refs. 30 and 31, and have recently been reanalyzed in the 
present context in Refs. 32 and 33. However, explicit model dependent threshold 
corrections for the SU(5) Yukawa couplings have not been computed and stud- 
ied. To properly generalize the analysis of Weinberg and Hall' 22 ' 231 to construct 
an effective MSSM, one must gauge-fix the superfield Lagrangian in a so-called 
supersymmetric S-covariant gauge so that the low energy gauge symmetry S is 
SU(3) x SU(2) x U(l). One essentially replaces derivatives by S-covariant deriva- 
tives in the gauge-fixing functional of the high energy theory. Details in the 
superfield formalism can be found in Refs. 34 and 35, but the corresponding ver- 
tices are not involved in the Yukawa threshold calculation. As mentioned above, at 
the high scale we can use the supergraph formalism' 361 to construct the effective ac- 
tion in Feynman gauge (£ = 1) obtained when the aforementioned heavy fields are 
integrated out. The one loop divergences are naturally regulated by dimensional 

[37] 

reduction (DR), which preserves supersymmetry at least up to one loop. 

[36 38] 

Due to the supersymmetric nonrenormalization theorems, ' the only modifi- 
cations to the parameters of the superpotential in the effective action arise through 
superfield wavefunction renormalizations. Hence if the original action has the form 




d 4 x / d 4 9(& a e 29V $ a + P[$}6 2 (9) + pi[$]6 2 (9) + ...), 



(3.4) 
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with $ a = being the decomposition into heavy/light superfields respec- 

tively, then the effective action (with an unconventional normalization of the ki- 
netic term), suppressing gauge kinetic and interaction terms, is 



S e -ff = / d 4 x / d A 6(Z ij <p\e 23V <P j + P'W 2 {6) + P' ] W 2 {6) + NR + . . .) . (3.5) 



As before, NR denotes effective nonrenormalizable interactions suppressed by in- 
verse powers of superheavy masses which, in the GUT case, will include operators 
responsible for proton decay. To evaluate Yukawa coupling corrections one need 
only consider supergraph two point functions with external light MSSM matter 
(s)fermion and Higgs doublet superfields and at least one superheavy superfield in 
the loop. Using the methods of Ref. 36 we find that the effective action contribu- 
tions to the light electron family, down quark family, up quark family and neutral 
Higgs superfield kinetic terms relevant for the Yukawa thresholds are 



+ EY^E C H^ + D C Y^DH^ ] + U C Y^UH {2) 
+ h.c. + NR+...) , 



where Z EiDtU 


are matrices in family and color space. Writing Z = 1 + K 


K E = 


- 3^ 2 A(o, Ml, 0)1 + 3(y (d) y W) T A(0, m 2 3 , 0) , 


K E c = 


- Qg 2 A(0, Ml, 0)1 + 3Y^Y^A(0, M\, 0) , 


K D = 


- 3g 2 A(0, Ml, 0)1 + (2y( u )ty(«) + Y^Y^)A(0, M\, 0) , 


K D c = 


- 2g 2 A(0, Ml, 0)1 + 2(Y (d) Y (d)t ) T A(0, M\, 0) , 


K v = 


- 3g 2 A(0, M 2 V , 0)1 + (2yMty(«) + Y^Y^)A(0, M\, 0) , 


K v c = 


- 2g 2 A{0, Ml, 0)1 + (y( u )ty(«) + 2Y^Y^)A(0, M 2 H , } , 0) , 


K H m = 


K H p = -3g 2 A(0,M 2 Hs ,Ml) + |A 2 | 2 (sA(0, M\, M 2 v ) 



(3. 





+ §A(0,0,M 2 ) + ^4(0,0,±M 2 )) , 
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where g is the GUT scale gauge coupling. Here the universal function A arising 
from the d dimensional loop integrals is 

/d d k 1 1 

where e = 2 — |. In the case p = we get 

A{^MlMl) = j^Q--F 1 {M%M 2 B )) , (3.9) 

where ^ = \ + ln47r — ^ E and F\ is defined in Appendix B. 

To determine the Yukawa threshold corrections we simply follow the generic 
component field analysis. First redefine the bare low scale effective light superfields 
so that they have a conventionally normalized kinetic term. This gives a relation 
between the bare high and low scale Yukawa matrices. The running relation is 
obtained by absorbing the divergent parts from the wavefunction renormaliza- 
tions using the definition of the high and low scale (3 functions. The resulting 
threshold matching functions relating the renormalized couplings in the region of 
the unification scale then depend only on the finite parts of the Zs. We obtain 
yM =y a (n)(l + A^ T ), where 



16vr 2 A^ = - y - [5Fi(M GUT , 0) + 3F l (M 2 Hs , M 2 GUT )) 
+ l(yi + f b )F 1 (M 2 H3 ,0) 

+ ^(3F l (M 2 H3 ,M 2 GUT ) + ^(MIO) + JL Fl( i M 2,0)) , 

16tt 2 A^ = - 9 -^{5F l (M 2 GUT ,0) + 3F l (M 2 H3 ,M 2 GUT )) 

+ (y 2 + ly 2 b ) Fl (M 2 H3 ,0) (3.10) 

+ ^{m{M 2 H3 ,M 2 GUT ) + §F!(M E 2 ,0) + ^F!(iM 2 ,0)) , 
16tt 2 A^ = - £ {m(M 2 GUT , 0) + 3F!(M 2 3 , M 2 GUT )) 
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+ l(y* + y 2 b )F 1 (M 2 H3 ,0) 

+ ^(3 Fl (M 2 H3 ,M 2 GUT ) + iF^MlO) + ^ Fl (lM 2 ,0)) , 

and where the (un)barred couplings are (low) GUT scale parameters (y& = y T ). 
The yb/Vr threshold matching function is therefore 

fy) = 1 + T ^(2, 2 (/,)(21n^ - 1) - 1^(2 In ^> _ 1)) , (3.11) 

and can be applied for /i ~ M GUT . 

We also quote the gauge matching conditions for minimal SUSY-SU(5) in the 
DR scheme [3 °' 311 : 

' ! -Ar T (/i), (3-12) 



where a G = g 2 /Arc and 



ArW = --ln^ + ^-ln M - 



3 



Af^) = - i In ^ + i In ^ , (3.13) 

Ar( , l )=-LA + l lD k + ± ln ^. 

3 VW 7T // 2tt /i 2vr /i 

In the interest of simplicity and minimality we do not include the effects of gravity- 
induced nonrenormalizable operators 1391 here. As noted in Ref. 33, differences of 
the Aj depend on M v and M E in the combination M GUT = (M 2 M S )? which we 
will take as GUT scale at which the matching functions are applied in the analysis 
to follow. 
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4. Electroweak and Supersymmetric Thresholds 



In addition to the GUT scale thresholds there are also thresholds at the elec- 
troweak scale and the effective supersymmetry scale M SUSY . At the electroweak 
scale, Mz, we integrate out the top quark, the weak gauge bosons and the SM 
Higgs. The definition of the inverse electromagnetic coupling ct(Mz) = 127.9 ±0.1 
and the Weinberg angle Sq(Mz) = 0.2324 ± 0.0003 includes one loop corrections 
from electroweak gauge bosons and the top quark 1401 for a pole mass Mto = 143 
GeV. This must be corrected for different top masses above Mz- The result is 1411 

1 3„ 2 /„. u/ 1 8 , M t \ 3 A to 2 p 
— ; r =-(1 - s 2 (M z )) —, r + — In — - ) -s— , 

ai(M z ) 5 V V Z!! \a{M z ) 9tt M tQ J ha{M z ) 

1 _ 2 / 1_ _8_ Mt_\ A$ P ( 41 ) 

a 2 {M z ) ~ S { Z '\^(Mz) + 97 11 M^J + l^Mz) ' 

1 1 , 1 , M t 

— In 



a 3 (M z ) a s (M z ) 3tt M z ' 

where Aj° p sa -0.92 x 10~ 7 GeV~ 2 (M t 2 - M 2 ) m and accounts for the quadratic 
top mass dependence of sin 2 6 W (M Z ) = s 2 (M z ) in Ref. 40. Below M z the "true" 
decoupled QED coupling is related to a(Mz) via 

1 ± ln ^»_i (1 + 211ll i^). (4.2) 



a-(M z ) a(M z ) 9vr M z Qtt M w ' 



We also incorporate the mass thresholds of Ref. 13 to determine the relation 
between the running fermion masses defined in the effective SU(3) x U(l) low energy 
theory and those above Mz- The effect of these thresholds is quite small given 
the smallness of a 2 {Mz) ~ 0.03 and the Yukawa couplings of the light fermions. 
However these results did not account for a heavy top quark. Essentially the only 
modification of the results is the matching function for the bottom quark mass. By 
properly integrating out the charged Nambu-Goldstone bosons in Feynman gauge 
one obtains a dependence on the SM top Yukawa coupling, yf M (yf M (fi) = rrii(fj,)/v 
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where v = 174.104 GeV and i runs over all SM fermions). The general result for 

[13] 

all massive fermions lighter than the top quark is 



where 



"»'""(/") = ™?"(f>)(l + A m") , (4-3) 



^ = |^((-IteW 2 + !(A) 2 )(ln^-|) + ^(ln^-I)) , (4.4) 



for all but the bottom quark mass and c = cosOw- Here g l VA = 2(g l L ± where 
9lr — T3LR~ s2 Q l an< ^ -^3 LR an< ^ Q l are ^ ne third component of weak isospin 
and the electric charge, respectively, for a given handedness of the i th fermion. For 
the different quark and lepton charge sectors one has 



9 V A = 


1, 


9v 


9A = 


-1, 


9 e v 


9A = 


1, 


9V 


9 d A = 


-1, 


9v 



- F. 



(4.5) 



For the case of the bottom mass one obtains 



*r = j^-^-hv + f )( ln ^ - i) + ^i 2 ^(<> ^?)) 



(y \^l%| 2 (^i(M^M, 2 )-lF 2 (M^M, 2 )) , 



(4.6) 



where the CKM matrix element \V t b\ ~ 1 and are defined in Appendix B. The 
term proportional to yf M is dominant. For pole masses Mt up to 200 GeV this 
effect corresponds to a shift of m^Mz) downwards by at most 0.6% from its value 
below Mz- 
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The pole masses for the heavier quarks, M q (q = c,b,t), are determined by 
finding a consistent solution to the three loop relation 



where 




with K c = 13.3 and Kf, = 12.4. In this analysis we will take the pole mass 
Mjy in the range from 4.7 to 5.3 GeV. Using the three loop running of the strong 
coupling, a s ] we find m b {M b ) = 4.06, 4.24, 4.43, 4.61 ±0.08 GeV and m h (M z ) = 
2.81, 2.96, 3.11, 3.26±0.16GeVforM b = 4.7, 4.9, 5.1, 5.3 GeV, respectively, and 
we have included the uncertainty for the strong coupling a s (Mz) = 0.120 ± 0.007. 
The pole mass for the r lepton is given experimentally by M T = 1.7771 ± 0.0005 
GeV,' 451 which, using the QED relation between the pole and running masses gives 
m T (M z ) = 1.7476 ± 0.0006 just below the electroweak threshold. With the loop 
effects of W and Z bosons included as in Eq. (4.4) one has m^(Mz) = 1.7494 ± 
0.0006 just above the threshold. 

The supersymmetric threshold is a potentially more important correction to 
to RG evolution of the couplings. It is sensitive to the details of the sparticle 
spectrum and in the case of the Yukawa coupling thresholds has quite different 
effects in different regions of parameter space. We shall first review the gauge 
thresholds and then give a simplified form for the Yukawa thresholds. The exact 
Yukawa coupling matching functions for the MSSM in the case of one light Higgs 
doublet are given in Appendix A. 

The exact form of the matching functions for the gauge couplings at M SUSY 
is well known. We will parametrize the gauge thresholds in terms of three mass 

1 We also include two loop QCD threshold effects' 44 ' which, however, are negligible in the 
mass range of interest. 
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rgj 

scales, Mi. The general result for the gauge threshold matching conditions is 

' ' -Af-»-Af«, (4.9) 



«i (/i) 

where af denotes the gauge couplings just above/below M SUSY . We convert from 



MS to DR couplings above M SUSY with the conversion factor given by Af R = 
-C 2 (Gi)/12n, where the quadratic Casimir C 2 (G) is N for SU(N) and for U(l) 
groups. The matching functions are 

Ar y (/,) = ^^ (p) ln^, (4.10) 

(p) 

where p runs over all sparticles integrated out near \i M susy . Here b\ is the 
contribution of sparticle p to the one loop coefficient b{ = 1, —3), (i = 1, 2, 3) 
of the (3 function for gi in the MSSM, 

n dgi bi 3 . . 

where t = In /x. The cumulative effect of these thresholds can be parametrized in 
terms of mass scales Mj via 

Ar» = i-(6,-6r)ln^, (4.12) 

where b? M = (|^, — ^P, — 7) are the corresponding coefficients in the SM and 6j — 
6f M = (|, ^,4). To get some idea of the dependence of the Mj on the sparticle 
spectrum we compute them for the case of separate degeneracies among squarks, 
sleptons, gauginos, higginos, and heavy Higgs particles: 

_9_ 11 J_ J_ 

Mi = m| B m| B m||- , 

_3_ _9_ _4_ _8_ J_ 

M 2 = mf 'mf 'mf r m,?~ r m'£ , (4-13) 

1 1 H W n 

1 1 
-^3 = rn~m~ . 

The complete Yukawa thresholds at M SUSY are quite complicated, however 
we can give some indication of their form and estimate their effects. The most 
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tractable form useful for making estimates occurs in the limit of small gaugino- 
Higgsino mixing (see Eqs. (A.30)-(A32)) and is given by 

Z/r(^) = yt(^)sin/3(l + A--), 

V§ M ( t M) = Vb (j*)<x* P(l + (4-14) 
y s T M M = yrM cos (3(1 + A s y -n , 

where 



167r 2 A^ sy ~ \yl[F 2 {ml T ,m\ ± ) + ^F 2 (to| ,m| ) + sin 2 f3 F 2 (m 2 H± ,0) 

+ ^ 2 2 ^ H ^ 2 tan/?^(G 2 (m|,m| ,m|,3)(O[ 1 ) 2 (4.15) 

- 2t 2 M 1 (fx H t a nP + A T )Y J G 2 (m\ , to? , m\ ) (i* ■ (r) ) 2 V 
16tt 2 A^ ~ |^ 2 ('-2m 5 ( / i H tan/5 + A)^G 2 (m|,m?,m|)(4.(6)) 2 

+y 2 ( -2/i„ tan/5 ^ (m t G 2 (m\ ± ,m\, m\)P ! id (t)P id (t) 

+ (ji H cot /3 + A t ) G 2 (m| ± , to? , to| ) (i* ■ (t ) ) 2 ) 

+ i(^F 2 (TO?,TO| ± )(0| 2 ) 2 + cos 2 /3F 2 (TO 2 f± ,TO 2 ))') (4.16) 

i 

+ l d fe ^K, to|j(o1 1 ) 2 + £ f 2 (to? , to| ) 

+ sin 2 /3F 2 (TO 2 ?± ,TO 2 )^ 



+i(7 2 2 U H 7W 2 tan (3 (^(m? , to| , to^ 3 ) (O^f 
+ G 2 (to|,to| ± ,to 2 ^)(01 1 ) 2 )) , 



20 



16tt 2 A^ ~ |y 3 2 f-2m^(m^ 2 (m|,m|,m|)Pi.(i)Py(i) 

hj 

+ {fi H cot (3 + At)G 2 (m|,m?,m| )(/*■(*)) 

/ * (4-17) 
+ 30? ( E F * (™| , m% ± ) (Qfx ) 2 + E ^ (m? , m% ) 

+ cos 2 /3F 2 (m^ ± ,m^) 
+ \vl (E ^2(m| , m| ± )(Of 2 ) 2 + sin 2 /? F 2 {m 2 H± , m 2 )\ , 

where t = tanOw, Hh is the supersymmetric Higgs mass parameter defined in 
Eq. (A.l) and A t ^ T and M.\, 2 are the soft supersymmetry breaking parameters 
defined in Eqs. (A. 2). All parameters as well as the function F 2 have an implicit 
dependence on the renormalization scale (~ M SUSY ) The function G 2 is positive 
and is defined in Appendix B. The P matrices involve products of the squark and 
slepton mixing matrices O t,b,T and are defined in Eqs. (A. 3) and (A. 22). Note that 
all terms proportional to the products PP' drop out in the limit of no squark or 
slepton mixing. To put these into perhaps more conventional notation in terms of 
the left-right squark and slepton mixing angles 6t, &b an d 6 T defined as in Eq. (A. 3), 
use the general relations 

G 2 (m 2 , m| , m\ ) P- 2 = \ sin 2 26 q (g 2 (m 2 , m| , m | ) + G 2 (m 2 , m| , m | ) ) 

i,3 

+ \ cos 2 26 q G 2 (m 2 ,m\,m\ 2 ) , (4.18) 
G 2 (m 2 , m\ , m| ) P[- P {j = \ sin 26 q (g 2 (m 2 , m| , m| ) - G 2 (m 2 , m| , m| ) ) , 



E 



where q — t,b, r. The explicit form of these angles in terms of the superpotential 
and soft supersymmetry breaking parameters of Eqs. (A. 1-2) can be found in Ref. 6, 
for example. 
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For small tan (3 one may further neglect squark and slepton mixing except for 
the stop mixing effects. In this case the gluino contributions to the quark wave- 
function renormalization and to the vertex diagrams of Fig. 1 give the dominant 
effect. Note that for ji H < one typically finds larger stop squark mixing, so that 
the threshold effect in yt should be largest in that case. 

At large tan/3, where the corrections to y& and y T tend to be more signifi- 
cant, the sbottom and stau mixing must be included. In this case enhancements 
proportional to tan (3 occur from vertex graphs containing gluinos, charginos and 
neutralinos. The gluino contribution from Fig. 1 dominates the correction to yj, 
for large tan (3. The chargino and neutralino contributions of Fig. 2 are also large 
in this regime. Depending on the relative signs of fi H and At, the part of the fi- 
nite diagram of Fig. 2 proportional to p, H Aty\ tan (3 can either enhance or decrease 
the gluino contribution. The divergent wino exchange graph tends to depress the 
gluino contribution slightly. The enhanced y T vertex corrections come from the 
order g^t&nft contributions of Fig. 3. In this case the bino exchange graph can 
be important due to the maximal hypercharge of the leptons. It is opposite in 
sign to wino and Higgsino vertex contributions. In analyzing the possible large 
tan f3 enhancements, it is important to note that the sign of most of the relevant 
vertex contributions is controlled by the sign of fi H (or ji H and At) and the var- 
ious large contributions can appear with opposite sign. In addition, if we take 
M SUSY — Mz then 1 the wavefunction renormalization contributions tend to be 
positive. In light of these facts a complete analysis is required. We will discuss the 
numerical significance of these thresholds in Section 7. 



1 The choice of M SUSY is not important when the full threshold corrections are included since 
they incorporate the effects of one loop running between different effective supersymmetry 
scales. 
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5. Semi- Analytic Unification Analysis 



We outline a semi-analytic analysis of gauge and Yukawa coupling unification. 
In particular, we discuss the solutions for M GUT , M fl3 , and m t (Mz) obtained when 
the P functions for the gauge and Yukawa coupings are integrated at the two and 
one loop level respectively. For the gauge couplings we must approximate the 
two loop Yukawa coupling effects and for the Yukawa couplings we can give only 
approximate solutions to the one loop P functions. We then look at the effect of the 
one loop threshold matching functions on these solutions. This of course can only 
give some idea of the dependencies of these masses on other parameters since a full 
analysis with one loop thresholds requires integration of the two loop P functions. 
The full numerical analysis will be given in Section 7. 

We use the condition of gauge coupling unification to determine part of the 
GUT scale spectrum from deviations from the naive unification condition. Thresh- 
old conditions at the GUT scale relate the mismatch of the gauge couplings to 
masses of superheavy particles. At the two loop level we shall consider the sce- 
nario, earlier considered in Ref. 33, in which one predicts M Hs , M GUT = (M 2 M S )3 
and the coupling constant at the GUT scale, a G , for fixed a, sin 2 6\y and a s at 
Mz- This scenario, which at one loop essentially reduces to determining the effec- 
tive scale M SUSY rather than M fl3 , is sensitive to the large a s error bars as well as 
boundary corrections at Mz and M SUSY , particularly since M H[i itself only appears 
in threshold matching functions. 

We also consider the implication of the Yukawa coupling thresholds derived in 
Sections 3 and 4 to the minimal SU(5) prediction yb(M GUT ) = y T (M GUT ). One may 
either fix M t to predict the pole mass Mf, (or the running mass m^Mz)) or use 
the accepted range of allowed M& values to predict M t (or m t (Mz)). In the former 
case the threshold corrections are straightforward to estimate and constitute a 
significant correction to which should be included into the analysis of Ref. 8 
(These Yukawa corrections essentially belong to their correction parameter p F .). 
In the latter case the analysis is more difficult and must be done numerically, 



23 



however we will see how both the gauge and Yukawa coupling thresholds filter into 
a semi-analytic expression for mt{Mz). 

In the first approximation M GUT , M Hs and a G are determined entirely by gauge 
coupling unification. They depend on Mt only through two loop effects in the gauge 
running and through threshold effects, primarily at the weak scale. The solutions 
to the two loop gauge coupling f3 functions generalizing Eq. (4.11) are 

1 1 , (bf M - bi) , M SUSY , bt , M GUT 



ai(M z ) a G 2tt M z 2n M z y ' 

where Aj = Af M +Af USY +Af UT and the 9i are the two loop contributions, including 
the effects of the top Yukawa coupling. By utilizing the threshold corrections given 
in Eqs. (3.12), (3.13), (4.1), (4.9) and (4.12), one obtains the solutions 

t G = In— = t + A tc , 

t H = ln ^ = ^ + A te , (5.2) 
1 1 



a G a G 



q + A QG , 



where the naive one loop results are 



f0 = 7T / l-2s 2 _ 2 1_\ _ 2 ^ M SUSY 



6 V d 3a s J 9 M z 

-1 + 6s 2 10 \ 5 M SL 

2 V ~~ol 3a7 y + 6 n ~M Z 

1 + 14s 2 1 79 M SUSY 
+ tt; 1- tttt- hi ■ 



£=£[-^^-^)+^^> (5-3) 



a° 12a 18a s 36tt M z 



and the threshold and two loop correction terms are 

= §b.^ + C -£<-«. + ",-2*) , < 5 - 4 ' 

7Q W 1 

A QG = — In --^ + <5 SM ! + « + 5 GU -T + —(-501 + 39^2 + 2# 3 ) , 
36tt M sc/Sy q g «c 36 V ; 

where s = sin9\y, 8°* = jg and 8 D -i = — jjg^- We will use superscript 0's to denote 
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the one loop results which use a naive step approximation to the thresholds. In 
the above all low scale parameters are implicitly evaluated at Mz- Note that the 
two loop contribution of the top Yukawa drops out of the solution for M GUT while 
if we estimate constant yt ~ 1 then this effect changes t H by of order -^yft G 0.4 
and \ja G by j^Vtta ~ 0.2, respectively. 

The weak scale threshold corrections are 

fiSM _ fitop;s 2 _|_ fitop;a _|_ ^top;a s 

= _^ + ± (1 - 2s2 ) ln ^i-l, n ^, 

3 a 27 V ; M t0 27 M z ' 

A " M t 5. M t (5-5) 
= 37r — - h -(— 1 + 6s J In — In 



a 9 V ' M m 9 M z ' 

( jSM i _ fitop^s 2 _|_ ^top;a _|_ ^top-^OLs 



+ — (-1 + 14s 2 )ln— ^ + — In. 



6 a 27vr y M t0 54tt 

The residual GUT threshold due to the undetermined value of M v is 

, „ = _5 ln «^ 5 

q g 7i M GUT 3tt M e v 7 

The three effective scales replacing the naive scale M SUSY are 

_ 25 25 /m~\ III 

£W = M 1 16 M 2 16 M 3 ~ (-i) mlm? , 

V SUSY = M x 4 M 2 4 M 3 " 4 ~ (-*) ( — —] m* m^ , (5.7) 

Vm g -/ V trig J h h 

25 325 _4_ /Uln\ 36 22 12 JL J. 

W S[/sy = M, 316 M 316 Mo 79 ~ (— M J m»»- 79 . 

1 z o \Tflq/ W H H y 

We can approximate a range for these effective scales in the context of the canonical 
soft breaking of supersymmetry induced from a hidden sector of N= 1 supergravity. 
In this case there is a common SU(5) invariant gaugino mass Mi, a common soft 
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super symmetry breaking scalar mass Mq and a common trilinear scalar coupling 
A (see Eq. (A. 2)) at M GUT . Then there is a simple one loop relation between the 
gaugino masses: 

a 3 (M sc/sy ) 

Assuming ~ rriq we obtain the approximate relations 
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Ususy ~ 0.53 wig , \/ 5lw « 0.08 m^m^ , W SUSY « m 2 ~m^m^ . (5.9) 

Alternatively, in the no-scale supergravity case (Mo = A = 0) one has ~ ~ 
3m^ ~ 3m^. We then have instead the approximate relations 

U 3U3 Y ~ 0.80 , V SUSY w 0.06 m|m| , WW « 0.70 mjm| . (5.10) 

As pointed out in Ref. 33, the supersymmetric threshold corrections to M GUT de- 
pend mainly on the wino and gluino masses while M Hs is most sensitive to the 
Higgsino mass. The corrections to a G depend mainly on M2 and are most sen- 
sitive to the squark masses. If we allow a range of 100 — 1000 GeV for m# and 
mg while the gluino mass is taken to range between 120 — 1000 GeV, 1 then for 
degenerate squarks and sleptons 64 < U SUSY & 530 GeV and 8 ^ V SUSY ^ 80 
GeV, and bear little relation to the naive Mz & M SUSY ^ 1 TeV. We also take 
100 ^ W SUS y ~ M2 ^ 700 GeV in this case. Note that the approximate upper 
limit on the sparticle masses of 1 TeV is the generic order of magnitude expectation 
if supersymmetry is to avoid introducing a new low scale hierarchy problem. In the 
no-scale case the corresponding limits are 96 ^ U SUSY & 800 GeV, 6 ^ V SUSY ^ 60 
GeV, and 73 < W SUSY ~ M 2 < 630 GeV. 



1 The lower limit on the gluino mass is near the experimental lower bound 1461 if the light 
gluino window is closed. 
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Before turning to the M t solution we give an estimate of the relative importance 
of the different corrections to t G and t H . The relevant low energy parameters 

[40,47] 

are 

M z = 91.187 ±0.007 GeV , 
a _1 (M z ) = 127.9 ±0.1 , 

(5.11) 

a s {M z ) = 0.120 ±0.007 , 
si = 0.2324 ±0.0003 , 

where, as in Ref. 41, we denote by so the sine of the Weinberg angle for a central 
value of Mto = 143 GeV with the quadratic top mass uncertainty removed and 
treated as a threshold correction as discussed in the previous section. The current 
best fit value of a s from LEP and all collider and neutrino experiments is a s = 
0.120 ± 0.006 ± 0.002 m where the second set of error bars corresponds to a SM 
Higgs mass in the range 60 < < 1000 GeV and we have included a conservative 
error in (5.11). 

For central values of the parameters the one loop solutions are t° G = 32.9 
corresponding to M GUT = 1.83 x 10 16 GeV and l/a° G = 24.5 for M SUSY = M z . We 
should note that the one loop solution, M® s given in Eq. (5.2), is only a definition 
and is related to the scenario is which the effective scale M SUSY is determined by 
the requirement of gauge coupling unification. We may write 



t HX = In — = In — jr- 2 - H In ■ 



M H3 , Mi 19, T, 

^GUT 

2tt 

+ y(50i - 12^ 2 + 7^ 3 ) + + . 



(5.12) 



M° 3 19 M SUSY 2tt/3(1-5s 2 ) 7 \ 19 M su£ , 



where 



In— f^ = — ln^^- = -— " v - J + — +— ln^^, (5.13) 
M GUT 18 M% SY 9 V a a s J 18 M z ' v ; 

and M® USY , the naive one loop value required for gauge unification, is 79.2, 7.3, 
0.9 GeV for a s {M z ) = 0.113, 0.120, 0.127, respectively. The scale T SUSY = 
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_ 25 100 _ 56 

M l 19 M 2 19 M 3 19 , using the notation of Ref. 48. For M SUSY = M z this gives 
M° 3 w 2.6 x 10 17 GeV for a s (M z ) = 0.120, however, M„ 3 can only be reliably 
determined at the two loop level. 

Next we quantify the correction terms to these naive estimates. The current 
experimental limit on the top quark mass from fits to electroweak data is Mt = 
164+^21 GeV,' 471 where the central value is for a Higgs mass of 300 GeV and the 
second set of error bars corresponds to 60 < < 1000 GeV. In the MSSM a 

■ -1 rj . £i [42] 

range 50 < < 150 GeV is more appropriate, giving Mt = 143_ 1 g_g GeV J The 
current lower limit from the D0 experiment is Mt > 131 GeV (95% C.L.).' 491 We 
will allow Mt to range from 130 to 200 GeV, and obtain the following ranges for 
the top dependent electroweak corrections: 



(-0.04, 0.24) 5l° G P ' a = (-0.008, 0.03) <5g p;Qs = (-0.013, -0.03) , 
(.39, -2.2) 5l° H P;a = (-0.02, 0.06) 5^ = (-0.20, -0.44) , (5.14) 
(0.05, -0.27) 5 top i a = (-0.005, 0.02) 5 top i as = (0.002, 0.005) , 

where the first(second) entries are for the lower(upper) limit for M t . The influences 
of corrections from thresholds, two loop terms and the experimental error bars in 
(5.11) are summarized in Table 1. The top dependent thresholds in (5.14) are 
summed together with DR conversion factors in the SM threshold entries for the 
above range of Mt values. The MSSM threshold entries include the supersymm- 
metric threshold as well as the effect of running SM parameters from Mz to M SUSY 
for the previously mentioned ranges taken for U SUSY , V SUSY and W SUSY . The upper 
range corresponds to the degenerate squark-slepton case while the lower entry is 
for the no-scale case. The a s error bar entries are given first for the lower, then 
the upper limit on a s . The percent deviations given are relative to the naive one 
loop predictions given above. We see that M H3 decreases with increasing M t and 
increases with increasing a s and Higgsino mass (V SUSY ). The effective GUT scale 
increases with increasing Mt and a s and decreases with increasing gluino mass 



<rtop;s _ 

°t G 

ctop;s 2 _ 

tH 



(Ususy)- The magnitude of the effect on M GUT due to each of these parameters is 
comparable, ranging up to 40%. 



Correction 


Parameter Variations 


% Deviations 




h H 


Sa G 1 


M GUT 


M »3 


i 


SM (top) threshold 


(+0.008, +0.35) 


(+.06,-3.0) 


(-0.03,-0.36) 


(0.8,41) 


(6,-95) 


(-0.1,-1.5) 


MSSM threshold 


(+0.08,-0.39) 
(-0.01,-0.48) 


f-2.0,-.H) 
(-2.3, -.35) 


(+0.07, +1.6) 
(-0.17, +1.5) 


( 8,-32) 
(-1,-38) 


(-86,-10) 
(-90,-30) 


( 0.3,6.5) 
(-0.7,6.1) 


Sq error bar 


±0.04 


±0.36 


±0.05 


±4 


(-30,43) 


±0.2 


a error bar 


±0.03 


±0.06 


±0.02 


±3 


±6 


±0.1 


a s error bar 


(-0.19, +0.17) 


(-2.8, +2.5) 


(+0.03, -0.03) 


(-17,18.5) 


(-94, +1120) 


±0.1 


two loop (y t = 0) 


0.18 


-3.8 


1.1 


20 


-98 


4.4 


two loop (y t only) 





0.4 


0.2 
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Table I. Corrections to M GUT , M H:j and 1/a, 



Note that the A 2 correction dominates the SM threshold effects and, as ex- 

s ' 

pected, the "correction" terms for M Hs can be large. The a s error bar correction 
is the largest, changing M H3 over three orders of magnitude. The MSSM threshold 
correction to M H3 was included in the analysis of Ref. 33, but the top dependent SM 
threshold was not discussed. These two effects can be comparable in magnitude, 
with the SM effects changing M„ 3 by up to two orders of magnitude, correspond- 
ing to a change in the proton decay lifetime of four orders of magnitude. Clearly, 
when we also solve for Mt from the condition of Yukawa unification, the solution 
for M Hg will be substantially correlated with the Mt solution. The tendency to 
predict large values of Mt means that the effect of the SM threshold will generally 
be to depress M H . V As we shall see in Section 6, the opposite is true for the lower 
bound on M Hs from proton decay. The combination of these results may be useful 
in strengthening proton decay constraints on the parameter space. [ ° 01 

We next turn to the Yukawa unification predictions for either Mt or m^iMz)- 
We discuss both cases with emphasis on the former. The one loop (3 functions 
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for the Yukawa couplings are not solvable analytically, however we shall give some 
semi-analytic solutions. The one loop Yukawa (3 functions for the third generation 
in the MSSM are 

dy a 



Py a = 



Va 



dt -^-Wi+tcfWp), 

where a = (t, b, r) and b a p and c a i are given by the matrices 



(5.15) 



/6 


1 






' 15 


3 


16\ 
3 1 


1 


6 


1 


and 


7 
15 


3 


16 
3 


\o 


3 


4; 




V £ 

\ 5 


3 


0/ 



(5.16) 



respectively. The analogous matrices in the SM, b^Z and c^f, are given by 



' 2 


3 
2 


1\ 




' 20 


9 

4 


8 \ 


3 
2 


9 
2 


1 


and 


1 

4 


9 
4 


8 


\3 


3 


2 / 




V 5 

\ 4 


9 
4 


0/ 



(5.17) 



respectively. At one loop, in the naive step approximation, the SM Yukawa 
couplings are matched to those in the MSSM at M SUSY via y t = yf M / sin (3, 
Ub — yf, M I cos P an d Vr = Ut M 7 cos Pi where tan (3 = ^ is the ratio of the VEVs 
of the two Higgs doublets H\ and if 2 in the MSSM (see Appendix A). Using 
Eqs. (5.15), the one loop solution for mf,/m T at Mz can be written as 



m b (M z ) _ o o o S M-§ osm§ Osm-I o 



(!)*■ 



(5.18) 



where 



A SM 
A b/T 



A 



a?(M 



SUSY J 



(A{M Z ) 

,,0/ 



12 ' dj(M gt , 5Y ) 
a 3 °(M z ) 



a5(M st/sy )/ V a 0(M s 



a. 



(5.19) 
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and 



fir 

B° a SM (M z , M SUSY ) = exp f -^-2 J dt (y° a SM ) 2 



B° a (M SUSY , M° GUT ) = exp ^- JL Jdt (y°) 2 j 



(5.20) 



Here the parameters t s = In M ^ Y , t z — 1 and i° is given in Eq. (5.3). On the 
other hand the full two loop solution including the Yukawa thresholds of Sections 3 
and 4 is 
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m b(M Z ) _ , SM . SM-l SM% SM-l ( B h 

- A b/r A b/T B t B b B T B t ^-jO b/T (52i) 
x(l + A^ + A^ + A^), 



where 



and 



= f Oi (M SUSY ) yir / a 3 (M SUSY ) yr 
h ' T \ a+(M z ) J \ a+(M z ) ) 

/ a{{M GUT ) \~I / 03 (M GUT ) \~| 
fe/r ~ L+(M S[7Sy )/ L+(M S[W )/ 



i?r(M+,M- sy ) = exp 



(5.22) 



B a {M+ USY ,M GUT ) = exp 




(5.23) 



The thresholds corrections A b / T are contained in Eqs. (3.11), (4.6) and (4.15) - 
(4.17). The two loop corrections are contained in Q b / T which we will take to be 
unity in this section. The ± superscripts indicate parameters just above/below the 
scale at which they are evaluated, i.e. on either side of the threshold boundary. We 
call the threshold corrections A b / T in (5.21) the direct threshold corrections. There 
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are also indirect threshold corrections which arise in rewriting the the quantities in 
(5.22) and (5.23) in terms of the naive one loop results. The discontinuities in the 
other parameters at the threshold boundaries indirectly filter into the result for 
m^/rrir. These indirect threshold effects will primarily involve the gauge coupling 
thresholds. 

For the indirect threshold effects in the A h j T parameters we find 



A 



b/r 



A 



b/r 



1 + A 



b/r 



A°b /r (l + ^/r) 



(5.24) 



where 



1 + K/r 



10 

l + a° 1 (M z )A s 1 M \ «/ l+a°(M z )A§ 



l + a\(M SUSY )^l M ) \l + al(M SUSY )kf 

l + qO(M s , sy )(Af* + Af^y 
l + ag(A aG -Af^(M GUT )) 
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X 



l + c^(M sc/sy )(Af M + A 



SE/SY + A 



3 > 



l + a°(A 



Af T (M GUT )) 



(5.25) 



Later we will also need 



A SM 



A° t SM {l + A? SM ) , 



(5.26) 



where 



A SM 



_ ( <Xi(M SUSY ) \ m. / a Q 2 {M SUSY ) yTs ( a%{M SUSY ) \- 7 
~ \ o?(M z ) > \ a\{M z ) J \ a°(M z ) ) ' 1 J 



22 ,„Q 



al{M SUSY )- 
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l + a?(M z )Af M \&/ l + a"(M z )Af 



l + a°(M sc/sy )Af 

1 + a°(M^)Af M 



X 



l + a°(M sc/sy )Af" 



l + a»(M SUSY )Af 

4 
" 7 



(5.28) 



for Af UT 



The threshold correction to a G is given in Eq. (5.4). In the expressions 
of Eq. (3.13) we shall insert the one loop expressions for In M H3 /M GUT given in 
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Eq. (5.13), consistent with the scenario considered here. Note that terms propor- 
tional to In M V /M GUT cancel in the differences A aa — Af UT so that these threshold 
corrections are computable given the sparticle spectrum. In all cases we find that 
A^ M and Af SM are negligible (< 0.1% for M SUSY up to 1 TeV) due to the absence 
of large logarithms. 

The indirect threshold corrections to the B parameters are primarily of second 
order compared to the direct thresholds in Eq. (5.21) and are more difficult to 
estimate analytically. However, one can account simply for the threshold effects in 
t G given in Eq. (5.4). We may write 

B **Bl(l-^A tG yty , (5.29) 

although the full corrections should be computed numerically. The threshold cor- 
rections to the B SM can be neglected. We can now rewrite Eq. (5.21) in terms of 
the naive one loop result 



m b(M z ) _ mg(M z ) Y 



(5.30) 



This may be used to determine m\,{Mz) using the Yukawa unification condition. 
However it is important to note that such a solution will be sensitive to potentially 
large threshold corrections both at M SUSY and M GUT) the former occurring for large 
tan /3 and the latter occuring for large allowed splittings M v ^> M E or M H . A ^> M GUT 
and large y t (see Sections 3 and 4). We shall discuss these corrections further in 
Section 7. 

We focus instead on using Eq. (5.21) to solve for mt(Mz). Here we use a good 
approximate analytic solution for y t for small to moderate tan/5 (see also Ref. 51). 

B t (M+ USY , M GUT ) ~ (1 - y?(M+ SY )K t (M+ USY , M GUT ))~& , (5.31) 
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where 



G 

K t (M+ SY ,M- UT ) = ^2 y dtA-\t s ,t) , (5.32) 



and 



A ,.+ f -S _ ( a i [ M GUt) \ 198 / a 2 [Mgut) \2/ a 3 [M GUT ) \ ~9 

t[s, G) y*t(M SUSY )) \4{M SUSY )) \ai{M SUSY )) ' ^ 
We then solve (5.21) for yt(Mz) using B t = k(%-f where k = fc°(l - 5k) and 



m T {M z y b ' T h / T> 1 b (5.34) 
Sh ~ A SM -i- A SUSY 4- A GC,T -i- A aSM 4- A a 



The solution for mt(Mz) is 

m t (Mz) = m° t (M z ) (l + A^ y + A? SM - ^A„ t 

+ C(5k + 3(A Bb -A BT ))^ , (5.35) 

where 
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m° t (M z ) = v sin (3A° t SM B° t ' SM > B° b SM *B° T SM ^ v ' 0/ , (5.36) 

f = 174.104 GeV and we have inserted the appropriate threshold corrections to yt- 
The factor C is given by 

c = T^v» ' (5 ' 37) 



1 - (A;°) ] 



Also, we have introduced a threshold correction factor A^ t , defined by K t 
K®(1 + Ajc t ), which must be determined numerically. 
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Next we evaluate some of the parameters appearing in Eq. (5.35). In particular 
we estimate the corrections arising from the gauge coupling thresholds. We will 
evaluate the Yukawa thresholds from supersymmetric particles of the MSSM and 
SUSY-SU(5) in the numerical section to follow. Table 2 gives values for the correc- 
tions to Vb/UT and Kt from gauge coupling thresholds which affect M GUT and a G . 
Results are given for tan/3 = 1.5 for explicit two loop Mt solutions. The results 
are insensitive to Mf and tan/3 (a two loop effect in the gauge coupling evolution). 
The gauge coupling sparticle spectrum parameters are taken to be M\ = 368 GeV, 
M 2 = 256 GeV and M 3 = 375 GeV. 



M SUSY 


M t 


A?/ 

b/T 


A* 


91.2 


157.4 


-0.06 


-0.03 


500 


173.9 


0.02 


0.003 



Table 2. Corrections to Mt from Gauge Coupling Thresholds in 
the MSSM (a s (M z ) = 0.120, M b = 4.9 GeV, tan/3 = 
1.5). 

The correction factors A# t are for central values K® = 0.77, 0.85 for M SUSY = 
Mz, 500 GeV, respectively. The correction term Ag b — A^ r in Eq. (5.35) is 
negligible except for large tan/3 in which case it has a value of —0.01 for M SUSY = 
Mz- For low tan/3 we also find m®(Mz) = (195.7 GeV) sin/3, consistent with the 
infrared quasi-fixed point solution' 52 ' 531 for mt- This behaviour is modified for larger 
tan/3 since the factor of B®/B® is proportional to exp(sec 2 /3). For low tan/3 the 
fixed point solution implies that this factor is very close to I. Thus the factor 
C multiplying the yj,/y r threshold corrections in Eq. (5.35) is determined by k° 
only. For the M SUSY = M z and tan/3 = 1.5 case in Table 2 we find k° = 0.72 
and C = 0.12. Hence we see explicitly how the fixed point solution inhibits the 
influence of potentially large GUT or MSSM threshold corrections to yb/y T on the 
M t solution. In this sense the solution for m t (Mz) tends to be more robust than 
the alternative rrib(Mz) solution entertained by some authors. Of course, away 
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from the fixed point B®/B® tends to be larger than 1, and the factor C need not 
be small. Still, for moderate values of tan/5, one finds C ^ 0.3, and to the extent 
that the approximation of Eq. (5.31) is still valid, the robustness of the M t solution 
is maintained except for extremely large tan/?. 

6. GUT Scale Constraints 

In order the extract any meaningful information concerning the parameter 
space we must inpose some constraints on the superheavy sector beyond that of 
unification. Proton decay is the strongest constraint on the Higgsino color triplet 
mass due to the dimension 5 operators which it induces! 541 A recent thorough 
analysis by the authors of Ref. 33 suggests that a conservative lower bound is 



M H3 > 5.3 x 10 15 GeV . (6.1) 



This corresponds to extreme limits on the superparticles, i.e. heavy squarks of order 
1 TeV in mass and a maximal ratio of squark to charged wino masses, low tan/3 
and a top mass around 100 GeV. The proton lifetime for the expected dominant 
decay mode p — > K + h>^ is given by' 331 



rip -> K+Uf,) = 6.9 x 10 31 yr 



0.003 GeV 3 0.67 sin 2(3 M„ 3 



TeV 

x 



Pn A s 1 + y tK 10 17 GeV 

2 



-1 



m w (G 2 (m 2 w , m?, m|) + C7 2 (m^, m|, m\)) 



(6.2) 



where fl n is the relevant nuclear matrix element and lies between 0.003 and 0.03 
GeV 3 . The parameter y tK gives the ratio of the third to the second family contribu- 
tions. Since it involves undetermined complex phases from CKM matrix elements 
one must allow for possible cancellations between different families. A conser- 
vative value of 1 1 + y tK \ 0.4 can be taken to get conservative limits on M H . 
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The quantity A s is a short-distance factor incorporating the anomalous dimen- 
sions of the important dimension 5 operators as well as the renormalization of the 
strange and charm masses from Mz to M GUT . The current experimental limit is 
r(p — > K + u ll ) > 1 x 10 32 yrj 551 implying that the squared factor in (6.2) must be 
^O(l). 

We would like to note that the limit (6.1) can be substantially strengthened 
when combined with sparticle spectrum predictions and the large M t solutions re- 
quired by bottom-tau Yukawa unification. Note that A s is enhanced by up to a 
factor of 3 when such large values of yt are included in the running charm Yukawa 
coupling, decreasing the predicted lifetime by nearly an order of magnitude. In 
particular, we have an additional motivation to compute M H3 rather than a s (Mz) 
as a gauge unification scenario. By determining M Hs for each sparticle spectrum 
determined by one's favorite set of theoretical criterion, e.g. universal soft SUSY 
breaking, radiative electroweak breaking, etc., one can immediately apply strength- 
ened and correlated limits from proton decay to this spectrum. The well-known 
figure of merit is typically the factor R = m~/(m^ ± x 1 TeV), where rriq is a first 
generation squark mass. This factor obtains a lower bound from proton decay, 
while typical sparticle spectra generated with the standard theoretical assump- 
tions require R ^ 25 (^ 4 in the no-scale case assuming rrig ^ 120 GeV) for 
sparticle masses below 1 TeV. Even more of the parameter space of the MSSM can 
be excluded by this procedure! 501 For the present analysis however, we will consider 
the conservative limit of Eq. (6.1). 

One may also bound the adjoint and color triplet Higgs masses from above by 
the requirement that the superpotential couplings Ai,2 be perturbative (< 2 V /F) 

[1 331 

below the Planck scale. ' Of course, this is purely theoretical prejudice as there 
may indeed be new physics between M GUT and M Planck ~ 1.2 x 10 19 GeV. We there- 
fore use the perturbativity constraint only to give an indication of a realistic GUT 
spectrum and not a strict requirement. We do not impose fine tuning constraints 
which require that these couplings are not too small' 11 as such small couplings in 
the superpotential are technically (though not aesthetically) natural due to the 
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nonrenormalization theorem.' 381 

In the one loop approximation the GUT scale renormalization group equations 
for the Yukawa-like couplings, taking into account only third generation effects, 
are 

ft, =^>Mf M + 3A| - 30 9 2 ) , 

ft. =I^lM§A5 + f A3 - I ,/ + 3,? + , (6.3) 

1 ,24- 2 2 96 ox 

1 ,24,, 9 9 84 9n 

^ =T^2^( y A 2 + 3rf + lOy? - -g 2 ) , 

where we have corrected an error in the (3 yt of Ref. 33. Using the minimal SU(5) 
relations 

M H3 A 2 M E 2Ai 



(6.4) 



M v g M v g 

we find the following results for one loop perturbativity up to M Planck /\f&n for 
M GUT = 2 x 10 16 GeV and a G = 1/24.5: 

M H3 < 1.9M V , M E < 3.8M V . (6.5) 

If we require perturbativity up to 10 17 GeV instead, then the factors becomes 2.7 
and 5.3, respectively. The limits are obtained for each Yukawa coupling by taking 
all other Yukawa couplings to zero in its (3 function. The same constraints applied 
to the ordinary Yukawa couplings gives yt(M GUT ) ^ 1.5(2.1) and Db{M GUT ) ^ 
1.4(2.0), where the values in parenthesis correspond to perturbativity up to 10 17 
GeV. However in the analysis below we choose to be somewhat inconsistent by 
requiring only that these couplings remain perturbative up to M GUT . 



; J ,8 



7. Full Renormalization Group Analysis 



A complete two loop numerical analysis is necessary to properly determine the 
influence of threshold effects on the solutions for M t) M„ 3 and M GUT . We will 
consider the effect of turning on and off various thresholds as well as the effect of 
varying various parameters appearing in the threshold formulae. 

From Table 1 we saw that M H . } decreases with mt and a and increases with 
s 2 , a s and V SUSY (rrijj) while M GUT increases with mt and a s and decreases with 
s 2 , a and U SUSY (m^ , rrig) . Just considering gauge unification for the moment, 
we can quantify these dependencies numerically, incorporating all the gauge cou- 
pling threshold corrections. We numerically integrate the relevant two loop (3 
functions (see Refs. 13, 56 and 57 and references therein) and use a modified glob- 
ally convergent nonlinear Newton method to find the solutions for M H3 and M GUT 
consistent with the mismatch of the gauge couplings at M GUT . The results are 
indicated in Figs. 4-7, and highlight the various dependencies mentioned above. 
In all the figures we indicate the dominant effect of the uncertainty in a s (Mz) as 
well as the influence of large and small tan/5 for central values of sin 9\y(Mz) and 
a(Mz)- In Fig. 4 we depict the logarithmic dependence of M GUT on the gluino 
mass for 120 GeV < m~ g < 1 TeV for fixed M t = 165.4 GeV. For fixed a s we 
see that the variation of M GUT over the entire range corresponds approximately 
to a 50% deviation. Similarly, Fig. 5 shows the variation of M GUT with Mt for 
130 GeV < m t [Mz) < 200 GeV for fixed m~ g = 400 GeV. Here the rise in M GUT 
is determined by the quadratic dependence of sin 9w on M t and gives a 20-30% 
deviation over the entire range of Mt- Note that the curves cut off at different 
values of M t corresponding to the nonperturbative limit of yt- 

The logarithmic variation of M Hz with V SUSY ~ .08m ^ is shown in Fig. 6 for 
100 GeV < m k < 1 TeV and M t = 165.4 GeV. Note that only in the case of 
large values of a s do we obtain consistency with even the most conservative proton 
decay bound, M Hz > 10 15 - 7 GeV. Here we get an order of magnitude change in 
M H3 over the range of V SUSY . The more complicated dependence of M H3 with Mt 
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for 130 GeV < m t (M z ) < 200 GeV is given in Fig. 7 for fixed V SUSY = 45 GeV. 
The initial fall of M Ha with M t is due to the electroweak threshold while its rise 
for large M t is due to the effect of large y t on the two loop gauge (3 functions. The 
variations in M Hz are on the order of 30%. 

Before imposing the condition of Yukawa unification, gauge unification con- 
strains the solutions for M GUT and M H . 3 . For M t between 130 and 200 GeV and 
0.5 ^ tan/3 < 60 we find the ranges 

1.0(0.96) x 10 16 <M GUT < 3.9(3.7) x 10 16 GeV , 

(7.1) 

1.7(1.2) x 10 13 <M H[i < 1.6(0.9) x 10 17 GeV , 

for the case of approximately degenerate squarks and sleptons, with the no-scale 
limits given in parentheses. These results are comparable with those of Ref. 33. 
Note that the lower limit on M GUT and the upper limit on M H[i will be sensitive 
to the restricted range of values of M t which result from Yukawa unification. The 
application of the most conservative proton decay bound from the previous section 
can be translated into a lower bound on the strong coupling constant. By tuning 
the remaining parameters so that M Ha is as large as possible we find 

a s {M z ) > 0.118(0.119) , (7.2) 

where again the limit in parenthesis corresponds to the no-scale case. 

To get some idea of the size of the GUT scale Yukawa coupling threshold cor- 
rections, we combine the renormalization group constraints with the perturbativity 
constraints of Section 6. The perturbativity constraints of Eq. (6.5) together with 
Eq. (7.1) imply 

M v >7.0(6.1) x 10 15 GeV , 
M E <0.9(1.2) x 10 17 GeV , 

where the two values correspond to perturbativity up to Mp lanck /y/8n = 2.4 x 10 18 
GeV and 10 17 GeV, respectively. If we further assume that no superheavy masses 
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lie above the relevant Planck mass scale, we find in addition 



^ £ 1.2 x 10 7 . (7.3) 

Although this may appear to be an extreme choice there is nothing to rule it 
out phenomenologically, in fact such a case was considered recently' 581 where the 

[59] 

fine-tuning of Ai is used to solve the Polonyi problem. What we regard as a fine- 
tuning of parameters may indeed be fixed by some not as yet not fully understood 
quantum gravitational model. In any case, supersymmetry insures that any tuning 
of superpotential parameters is at least technically natural as they are multiplica- 
tively renormalized. From Eqs. (3.11), (5.34) and (5.35) we see that such a large 
splitting tends to increase the top mass prediction. In many cases no perturbative 
solution can be found, however, away from the infrared quasi-fixed point of yt, siz- 
able enhancements of the predictions can occur. In Fig. 8 we plot the Mt solution 
vs. tan /3 for various values of the ratio M v /M^ for M5 = 5.2 GeV, a s (Mz) = 0.120 
and M SUSY = Mz- The largest deviation from the naive case is a 15% increase 
in Mt and corresponds to A^_ T ~ 0.1. The sparticle spectrum parameters are 
fixed at rrig = M3 = 1 TeV and V SUSY = 80 GeV. The typical GUT masses are 
then M GUT = I0 16 - 1±ai GeV and M Hs = 10 15 - 7±a4 GeV for these solutions. We 
also plot the bound from the nonperturbative limit on y t , indicating that for this 
particular example, one is not in the domain of attraction of the fixed point. 

The GUT scale matching function for yb/y T of Eq. (3.11) can also be negative 
for large yt and M Hr We set M V /M E = 0.3, near the minimum required by the 
perturbativity argument. In Fig. 9 we give an example for this case for M& = 
4.9 GeV, a s (Mz) = 0.127 and M SUSY = Mz- Using the same set of sparticle 
spectrum parameters as in the previous figure, we obtain GUT masses in the 
ranges M GUT = I0 16 - 2±ai GeV and M„ 3 = 10 16 - 8±0 - 3 GeV. The deviations here are 
smaller: A^_ T ~ —0.02 since M H:j cannot be larger than M GUT by more than an 
order of magnitude. The important thing to note is that the lack of knowledge of 
the GUT scale spectrum does not preclude extracting reliable predictions for Mt 
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over much of the parameter space and the uncertainties in these predictions can be 
reasonably estimated. In regions where predictions are not robust, i.e. away from 
the attraction of the infrared fixed point for y t , one must worry about 10 - 20% 
effects in the minimal model. 

For completeness we also depict the most significant uncertainty in the Mt 
solution: the error bars on a s and Fig. 10 shows the M t vs. tan (3 contours 
for different a s for the generic case M v = M E , = 4.9 GeV and sparticle 
spectrum parameters rrig = 400 GeV, V SUSY = 35 GeV and M3 = 375 GeV. In 
Fig. 11 we show these contours for different ranging from 4.7 to 5.3 GeV for 
a s (Mz) = 0.120 and the same remaining inputs as in Fig. 10. Note that for small 
Mb the curves are cut off due to the lack of a perturbative solution for M t for 
intermediate values of tan f3. Clearly, in the absence of additional theoretical cuts 
such as proton decay or radiative electroweak symmetry breaking there remains a 
large parameter space. Even if one imposes the restricted Mt limits from fits to 
electroweak data, one can find solutions for all regions of tan (3, at least for larger 
Mfr. A better determination of M& could eliminate the intermediate tan (5 solutions, 
although this subject is frought with theoretical uncertainties. 

Next we turn to a discussion of the Yukawa threshold corrections in the MSSM. 
Using the approximate matching functions of Eq. (4.15) - (4.17) we incorporate M t 
into the numerical routine together with M GUT and M Hs and compute the threshold 
effects on the top mass solution for representative sparticle spectra consistent with 
universal soft supersymmetry breaking parameters at the GUT scale as well as 
radiative electroweak symmetry breaking.' 61 We also consider the low, intermediate 
and high tan (3 cases separately. Table 3 gives the effect of the threshold corrections 
for low tan/9 = 1.47 and some sample soft breaking parameters including the 
generic no-scale case. We also consider both signs of \x H and take A = throughout. 
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Mo 


Mi 




A sus Y 


A SUSY 


A SUSY 
Vr 


Mt (naive) 


M t 


100 


140 


+375 


0.02 


0.006 


-0.006 


164.1 


167.6 


100 


140 


-375 


0.03 


0.03 


-0.01 


164.1 


170.5 





150 


+353 


0.02 


0.007 


-0.007 


164.2 


167.8 





200 


-466 


0.07 


0.04 


-0.01 


163.0 


175.4 



Table 3. Low tan/3 Threshold Corrections to Mt in the MSSM 
(tan/5 = 1.47, a s (M z ) = 0.120, M b = 4.8 GeV, A = ) 

All masses in the tables are in GeV. Note that the increased stop squark splitting 
in the case of fi H < increases the yt corrections. The SUSY thresholds increase 
the top mass prediction relative to the naive result which includes only the gauge 
coupling thresholds. The dominant effect arises from the gluino-induced vertex 
corrections of Fig. 1 (for fi H < 0) and the gluino contribution to the top quark 
wavefunction renormalization. These are also the dominant contributions to the y b 
threshold, although for /i H > the gluino-induced vertex and wavefunction renor- 
malization contributions almost cancel, while they add constructively for \x H < 0. 
Since yt is near the fixed point in this region, the order 5% corrections to yb/yr 
translate into a negligible effect on the M t solution as discussed in Section 5. 

For intermediate and high tan/3 the bottom and r Yukawa corrections be- 
come more important. For intermediate tan/3 = 15 and Mi = 100 GeV, Table 4 
summarizes the typical threshold effects. 



M 




A SUSY 

yt 


ASUSY 

Vb 


A SUSY 
Vr 


M t (naive) 


M t 


100 


164 


0.01 


-0.1 


0.025 


194.5 


178.3 


200 


316 


0.03 


-0.03 


0.025 


190.3 


187.5 


300 


462 


0.035 


-0.01 


0.02 


187.3 


188.7 



Table 4. Intermediate tan (3 Threshold Corrections to Mt in the 
MSSM (tan/3 = 15, a s (M z ) = 0.118, M b = 4.9 GeV, 
Mi = 100 GeV, A = ) 
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It is interesting to see the source of the enhanced y^ and y T corrections. For 
Mq = 100 GeV the gluino induced vertex of Fig. 1 contributes —0.13, while the 
chargino induced vertices of Fig. 2 with a stop squark internal lines contribute 
0.0175. Accounting for other small corrections leaves a 10% effect. On the other 
hand for Mq = 300 GeV and correspondingly larger squark and slepton masses, 
many contributions come into play. Now the chargino contributions add to 0.03, 
the gluino vertex correction is —0.09 and the total wavefunction renormalization 
contribution is 0.03. Accounting for heavy Higgs effects and the neutralino contri- 
butions leaves a net 1% effect in Clearly a complete analysis is important even 
for intermediate values of tan/9. The dominant corrections to y T come from the 
neutral and charged wino induced graphs of Fig. 3. 

For larger tan (3 = 40 and Mi = 100 GeV, we give typical threshold effects in 
Table 5. 



Mi 




ASUSY 

yt 


A SUSY 

Vb 


A SUSY 


Mi (naive) 


M t 


200 


205 


0.04 


-0.035 


0.04 


177.6 


154.5 


400 


450 


0.05 


-0.045 


0.06 


169.0 


78.2 


600 


660 


0.055 


-0.03 


0.065 


162.3 


45.2 



Table 5. Large tan/3 Threshold Corrections to M t in the MSSM 
(tan/5 = 40, a s (M z ) = 0.118, M b = 4.9 GeV, M = 400 
GeV, A = ) 

Here the heavier spectrum chosen enhances all the wavefunction renormalization 
contributions and the gluino vertex contributions to yj, approach 20%. For the 
positive fi H chosen however, this is partially cancelled as described in the inter- 
mediate tan (3 case. Even so, the effect on M t is dramatic. Here one is far from 
the attraction of the infrared fixed point. The typically large M t predictions in 
SUSY-GUTs normally arise because the QCD running of y b typically places it 
below y T at M GUT . A large top Yukawa coupling counters this behaviour. How- 
ever, here the large, order 10%, threshold correction increases y^ relative to y T at 
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M GUT and smaller values of M t are required. In fact, for some spectra the situa- 
tion becomes unstable in the sense that no positive yt can be found such that the 
Ub(M GUT ) = y T (M GUT ) unification occurs. 

Finally, we use the results of Section 5 to estimate the contribution of the 
gauge coupling thresholds to the Mt solution which are already incorporated in the 
"naive" solution in the above tables. Using Eq. (5.24) we find the ranges A^ r 
(-0.014, -0.03), (-0.03, -0.045), (-.044, -0.084) for the data of Tables 3, 4 and 
5, respectively. In the intermediate and high tan/3 regions this effect supplements 
the other potentially large negative Yukawa corrections to yb/Vr- There is also 
a threshold correction to Kt of Eq. (5.32) given by A#- ( pa —0.02 for the data of 
Tables 3 and 4 and in the range ~ (—0.025, —0.046) for Table 5. The corresponding 
fractional correction to Mt is — \k-K t an d enhances the generally positive direct 
Yukawa threshold correction to yt- The generic features of the analytic solutions 
are born out, at least in the low tan (5 region where they are valid. In particular, the 
gauge and Yukawa corrections to yb/y T both at M SUSY and M GUT do not directly 
feed into the M t solution. They are multiplied by a factor, C of Eq. (5.37), which for 
the data of Table 3 is C ~ 0.1. The fact that the corrections to M t are dominated 
by the direct threshold correction Ay USY for low tan (3 confirms this expectation. 



8. Conclusions 



We have given a complete treatment of Yukawa coupling threshold corrections 
in the MSSM and minimal SUSY-SU(5). We have applied these as well as the 
gauge coupling thresholds to both approximate one loop analytic and consistent 
two loop numerical solutions for the superheavy masses M GUT and M H3 and the top 
quark mass, M t in the context of gauge and Yukawa unification in SUSY-SU(5). 
We have highlighted the sensitivities of these solutions to the low scale parameters. 
The effective GUT scale was the most constrained and variations of a s (Mz), Mt 
and the gluino mass led to of order 40% deviations in M GUT . 
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We have also found limits on the colored Higgs triplet superfield mass, M H . V 
from an RG analysis of unification. This approach avoids ad hoc assumptions 
about the degeneracies of the superheavy particles. We displayed the sensitivities 
of M H[i on a s (Mz), M t , tan/5 and the Higgsino masses which are much larger 
than for M GUT . Using these results we found that a conservative analysis of the 
proton decay bound could be translated into a lower bound on the strong coupling, 
a s [Mz) *S 0.118 for any sparticle spectrum with masses less than order 1 TeV. We 
also concluded that M Hz must be lighter than 2 x 10 17 GeV. This in turn limits 
the size of GUT scale corrections to yb/y T so that these effects cannot produce too 
large a decrease in the solution for M t . We found also that away from the attractive 
region of the infrared quasi-fixed point of the top Yukawa coupling, there can be 
significant positive corrections to M t if there is a substantial allowed splitting of 
the superheavy vector and adjoint scalar masses, M v 3> M E . This typically occurs 
for larger values of M&, while for smaller values no perturbative solution can be 
found. For low tan/5 these corrections to yb/Vr have different effects on the two 
alternative solutions, Mt for fixed or m^Mz) for fixed Mt. In the former 
case the corrections are suppressed, while in the latter they feed directly into the 
uncertainty for m^(M^). 

We have given a detailed calculation of the Yukawa corrections in the MSSM 
and analyzed their effect on the Mt solution for various regions of tan (5 and rep- 
resentative sparticle spectra. We have shown, both analytically and numerically, 
the robustness of this solution for low tan /3, where the dominant corrections come 
from gluino-induced effects. For larger tan/5 we discussed the importance of the 
complete calculation to determining the potentially large corrections to yi, and y T . 
These corrections were particularly sensitive to the signs of fi H and Af. For all 
tan (3 we described the enhancing effect of stop squark mixing on Mt for both signs 

of fX H . 

We emphasized the effect of the a s and error bars on M t , by far the 
dominant uncertainty. With strong bounds from proton decay in the context of 
canonical supergravity induced soft supersymmetry breaking, this parameter space 



4G 



can be cut substantially. A critical analysis of the theoretical uncertainty in 
would also help reduce the parameter space, particularly the intermediate tan/3 
region, which is still allowed if one accepts larger values of M b ~ 5.2 GeV. In- 
termediate values of tan/3 must also be reconsidered when the supersymmetric 
Yukawa thresholds are included as these can reduce M% by 10%. We will discuss 
these issues and perform a more complete analysis of the MSSM and GUT Yukawa 
thresholds under more specific assumptions for the sparticle spectrum in a future 
work.' 501 
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APPENDIX A: Yukawa Threshold Corrections in the MSSM 



Below we give the complete supersymmetry breaking threshold corrections to 
third generation Yukawa couplings in the MSSM in the case in which intergenera- 
tional mixing is negligible in the squark and slepton mass matrices. The relevant 
factors as in Eq. (2.10) are the finite parts of the one loop wavefunction renormal- 
izations for the third generation fermions and the light SM Higgs and the Yukawa 
vertex corrections involving heavy sparticles. We give results for the scenario in 
which there is only one light Higgs below the scale at which the threshold conditions 
are applied. 

First we give some notation and conventions. Our conventions closely follow 
those of Refs. 60 and 61 which the reader should consult for more details. We 
define our superpotential in the MSSM as 

P = e tJ (Y d H[QW + Y e H\VE + Y U H 3 2 Q % U + ^ H H\H 3 2 ) , (A.l) 

where the hats denote superfields, ei2 = — 621 = 1 and Y u ^ e are the conventional 
Yukawa matrices, with yt, y& and y T the respective diagonal elements for the third 
family. The superfields Q, U, D, L and E have the hypercharge assignments of 
the corresponding quarks and leptons. The Higgs superfields H\p have hyper- 
charge y = —1,-1-1 and their scalar components acquire vacuum expectation values 
(Hi) = vi/V2 and (H 2 ) = V2/V2, respectively. A general set of soft supersym- 
metry breaking parameters is introduced without explicit regard to their origin. 
These include soft trilinear scalar couplings mimicking those of the superpotential 
as well as explicit mass terms for the gauginos, squarks, sleptons and Higgs fields. 

V sott = - e l3 {Y d A d H\QW + Y e A e H\V E + Y U A U H J 2 Q*U + n H BH\H 3 2 ) , 
Cff = \Mi\ b \b + \M^ W X W + iM^Xf (A.2) 

- mgQ^Q* - m 2 jjU*U - m 2 D D*D 

- mlL^V - m 2 E E*E - m 2 Hi H[*H{ - m 2 H2 H l 2 *H l 2 , 
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where we denote superpartners of ordinary particles with a tilde, and A is used to 
denote the Majorana gaugino fields. We also define squark and slepton fields with 
the same charge conventions as their partners: Q = (ul^l), U = u* R , D = d* R , 
Q = (j>L, &l) and E = e* R . 

The threshold corrections will in general involve squark and slepton mixing 
matrices for the third generation. For example, denoting the stop squark mass 
eigentates as tip, we can relate them to weak eigenstates {l,r by an orthogonal 
matrix O* such that O t M'~O tT = M~ where M~ is diagonal and 

h\_f cos9 t sm9 t \ (~t L \ _ ( t L 
h ) ~ \-sm6t cose t ) \t R J ~ \t R 

We define O b and O r similarly. The neutralino mass matrix is diagonalized by the 
unitary matrix Z such that Z*M x oZ~ 1 = M x o and relates the weak bino, wino 
and Higgsino eigenstates to the mass eigenstate Majorana fields x° via 

-i - 2 o (A. l) 

H\ = Z* 3 Xi , H 2 = Z* 4 Xi • 

The upper index on the Higgsino fileds is an S\J(2) L index. Finally we define the 
unitary mixing matrices for the charginos by U* M X ±V^ = M x ± where the charged 
Higgsino and wino are related to the Dirac mass eigenstates Xdi = (xfi ^-i) via 

Hl = V* 2 xt , Hl = U* 2X - , -iX+=V^xt, -tX w = U* lX ;. (A.5) 

Explicit forms for the mass and mixing matrices can be found in Refs. 6, 61, and 
62. 

We first give the finite parts of the wavefunction renormalization of the top 
and bottom quark and r lepton due to loops involving squarks, sleptons, gluinos, 
neutralinos, charginos, and the heavy Higgs doublet involving the charged Higgs, 
H^, the pseudoscalar A and the heavy scalar H°. We sum the coefficients of $Pl,r 




49 



in the effective action obtained by integrating out these fields and find (using the 
notation of Section 2) for the top quark, 



K t (L+R)= -^2^IX! F2 Kv m I) 

+y 2 (j2 F 2 (ml,m 2 4 )(O b l2 ) 2 \U j2 \ 2 + F 2 (m 2 H± ,m 2 ) sin 2 /?' 

i,3 hj 

+ cos 2 /3(F 2 (m|±, mj) + \{F 2 (m 2 H0 ,m 2 ) - F 2 (m 2 A , m 2 ))) 
+9l fe F 2 (m| , m 2 ± ) (O b n ) 2 \U 3l \ 2 ( A.6) 

+ \ J2F2(ml,mlo)((O t il ) 2 \Z j2 + I^) 2 



+ ¥(of 2 ) 2 i^ii 2 t 2 )) 

f2 ^ F 2 (m? , m 2 )Of 1 Of 2 Re(C/ i iC/; 2 ) 



+^fv / 2^F 2 (m|,m 2 o)Oi 1 0| 2 Re((Z j2 - %*)Z* 4 ) 



and for the bottom quark, 
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mf) 



Wt{^F 2 {ml ) m 2 ^){0\ 2 ) 2 \V j2 \ 2 + F 2 (m 2 H± ,m 2 ) cos 2 (3 
+y 2 (J2^2 (m\ , m 2 xf ) (Oj, ) 2 1 U j2 | 2 + ^F 2 (m? , m 2 „ ) | Z j3 \ 
+ sm 2 p(F 2 (m 2 H± ,m 2 t ) + \(F 2 (rn 2 H0 ,m 2 h ) - F 2 {m 2 A , m 2 h ))) 

+9 2 2 (J2 F2 K > m lf ) ) 2 1 ^ 1 2 (A.7) 
+ ^F 2 K i ,m 2 o)((O l b 1 ) 2 |Z j2 - I^tl 2 



+ |(Qfe) 2 l^il 2 * 2 )) 

^2 F 2 (m|, m 2 ± )Ol 1 0| 2 Re(^iy/ 2 ) 



where F 2 is defined in Appendix B and t = tanOw- The result for the r lepton is 

K(L+R) = -^ 2 [y 2 r(j2 F ^ m h m \f ) I Uj2 1 2 + E F 2 K » ^ ) I ^3 1 2 

V j i,j 

+ sm 2 p(F 2 (m 2 H± ,0) + \(F 2 (m 2 H0 ,m 2 T ) - F 2 {m 2 A , m 2 ))) 
+^fE F 2«^ 2 ± )l^l| 2 



+ ^^K.m^^) 2 !^ + Z^| 2 (A.8) 
+ 4(0[ 2 ) 2 |Z jl |V 
-^(V^F^m? ,m 2 )O[ 1 O[ 2 Re((Z j2 - %i)Z* 3 ) 
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where K a{L+R) = K aL + K aR . 

The neutral Higgs wavefunction renormalization contributions arise from char- 
gino and neutralino loop diagrams in general giving matrix contributions so that, 

Z H = . > ( A - 9 ) 



where 



1 + K ho 



V i,j 

+ 2\Aij\ 2 F 3 (m y ±,m y ±))] , 



hJ 



+ 2| J B ii | 2 F 3 (m Y± ,m y± )) ] , 



2 



h3 



167T 2 

+ 2Re(A ij B* j )F 3 (m x ±,m x ±)) J , 
and where 

A"j = — Q'lj cos/5 + S"j sin/3 , Aij = —Qij cos/3 — Sij sin/3 , 
B'/j = Q"j sin/3 + S'/j cos/3 , Bij = —Qij sin/3 + Sij cos/3 , 

with 1601 



(A.10) 



(A.ll) 



1 1 

Qij = ^VnUj 2 , Sij = -j=V i2 U jl , 

Q'^i[Z i3 (Zj 2 -Z jl t) + (i^j)}e i , 
S'lj = \[Z iA {Z j2 - Zj X t) + (i <-> j)]ei . 

The factor e« = ±1 is introduced' 631 to insure positive mass neutralino eigenstates. 
When Z is defined without regard to the positivity of the mass eigenvalues, the 
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elements of the matrix M x o are e^m^o (i = 1, ...,4). The function F3 is given 
in Appendix B. In general the matrix Zh must be diagonalized, Zjj = OhZnOj^, 
giving an effective light Higgs field: 

K ff = (Z d H )l((O h ) 2l H + (O h ) 22 h ) . (A.12) 

This results in a threshold correction to the neutral Higgs mixing angle a of Ref. 60. 
However, in the limit of small gaugino-Higgsino mixing considered later, these off- 
diagonal contributions vanish. 

We next give the one particle irreducible Yukawa vertex corrections. These 
contributions involve more complicated functions of the squark, slepton, neutralino 
and chargino mixing angles although many of them simplify in certain approxima- 
tions for the sparticle mass matrices. The greatest simplifications occur in the limit 
of no squark and slepton mixing and vanishing mixing between the bino, neutral 
wino and the Higgsinos. Except for the case of large tan /3 one would expect small 
L-R mixing in the sbottom and stau mass matrices. However due to the large top 
Yukawa coupling the mixing can be substantial. Also the lack of mixing in the 
neutralino sector is expected if the bino and neutral wino masses are dominated by 
their soft-breaking masses and if the Higgsino masses are determined by the pa- 
rameter ji H in the superpotential. This approximation is certainly appropriate to 
the case of radiatively induced electroweak symmetry breaking near the low tan (3 
fixed point 16 ' 641 and works reasonably well for larger tan/5. We quote the entire 
results as well as the simplified results in the limit of small gaugino-Higgsino mix- 
ing. Note these terms can give rise to substantial enhancements of the threshold 
corrections proportional to ratios of sparticle masses to m& or m T if the mixing 
and mass-splitting between neutralinos or same generation squarks and sleptons is 
large. 

The finite Yukawa corrections are denoted = Sy a /y a as in Section 2. We 
separately give the results for Higgs, neutralino, chargino and gluino corrections. 
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For the HqU vertex we obtain 

KJ = ff?(Higgs) + Kj(x°) + KJ{ X ±) + Kj{g) , (A.13) 
where the separate contributions are 

^(Higgs) = ^ (iyt^iimlm^) - F^rrfo 

+ Ml (-(cos 2 2(3 - 2 sin 2 2f3)G 2 (m 2 Ho , m 2 Ho , m 2 ) 
+ 12sin 2 ( 5cos2 / 9G , 2(m 2 7o ,m| o ,m 2 ) 
+ cos 2 2(3G 2 (m 2 A ,m 2 A ,m 2 ) 

+ 4 cos 2/3 sin 2 /3 G 2 (m 2 4 , m%, m 2 ))^ ( A ' 14 ) 

+y?^Fi(mf,m| r± ) 

+ Ml ((2c 2 - cos 2 2(3)G 2 (m 2 H± ,m 2 H± ,m 2 b ) 

+ 2 cos 2 2f3G 2 (m 2 H±) m 2 V) ml)) 



+ m x o m x o G 2 (mf. , m 2 o , m 2 o ) JVgjj (t)^ (A. 1 5) 

+ 2m t ^2 m x oG 2 (m 2 x o,m 2 .,m 2 k )N ijk (t) J , 

i,j,k J 
V i,j,k V 

+ \ ± ^^Kv ra ^y^w) ( A - 16 ) 

+ 2m b ^m x ±G 2 (m 2 x ±,m^,m 2 b JC ijk (t) J , 
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KUg) = — ^( -yff3^^X^G 2 (m|,m?,m|)G^(*) J , (A.17) 

and where the functions F{ and G{ are defined in Appendix B. The quantities 
Nfj^(t), N^ k (t) and N^^t) in the neutralino contribution are 



where 



and 



Kjk(t) = ytN} jk {t) + g 2 N? jk (t) + 92(f)Nf jk (t) , 

N t f k (t) = y t N} jk {t) + g 2 N? jk (t) + g 2 { 9 j)Nf jk {t) , 

N ijk (t) = yfN] jk (t) + ytg2~N 2 i jk (t) + gl~N% k {t) 
+ alifwUt) + gl{f?N% k {t) , 

Vt Vt 



N} jk (t) = Re[TlUt)T° Rij (t)(-Q'%cotp + S'%)} , 
Nf jk (t) = - Re[(R L * k W° Rij (t) + T° L * k (t)R Rij (t)) 

x(-Q'%cotp + S"%)}, 
Nf jk (t) = Re[ J R^(t) J R^.(t)(-g // * fc cot/3 + ^)] , 



(A.18) 



(A.19) 



N] jk (t) = Re[TlUt)T Rji (t)(Pjk(t) + Pjk(t)( - COt f + > 

ttii 

N 2 ijk (t) = - ReftWW + TlIMRriM) 

x(p Jfc (t) + p; fc (t) ( ^ C0 ^ + ^ ))], 

= ^(cot 2 /? - ^RefT^^t)^,^)^^)] (A.20) 

+ Re[^(^)^(t)(^(t) + ^(^)( ^ C °lf + ^ ))] , 

fill 



Nt jk (t)= -^(cot 2 /?-l)Re[(i?°yt)^^ > 



The quantities Nij k (t) are obtained by interchanging the L and i? subscripts in 
(A.19). The matrices R and T involve products of squark and neutralino mixing 
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matrices while the P matrices involve products of squark mixing matrices only. 
They are defined by 



R° Lij (f) = v^C&CTa, - e s )tZ] x - T 3f Z* 2 ) , 
RRijif) = V20{ 2 e f tZ jX t 

T Lij(f) = °i2 Z jQ ' 

T° Rij (f) = O^Z^ej , 



(A.21) 



and 



W) = o! 1 o f jl + o f a o f j2 , 

P! j (f) = h(0{ 1 f j2 + 0{ 2 0j 1 ), (A.22) 
Pij(f) = O&j^Tv - e f s 2 ) + 0{ 2 O f j2 e f s 2 , 

where / = (t,b,r) with charge ej and weak isospin T3J. In T® R the ( 3 ) indices 
correspond to f — ( & T ) . For the sneutrino case there is no z> r # so one must make 
the replacements P Vt = 1, P' Vt = and P" t — \. In all these quantities c = cos 9w 
and s and t are the sine and tangent, respectively. 

The quantities C^{t), C^ L k (t) and Cij k (t) in the chargino contribution are 



where 



C$(t) = y b Cl jk (t) + g 2 C? jk (t) , 
C$(f) = y b Cl jk (t) + g 2 Cf jk (t) , 
Cijkif) = y 2 c\ jk {t) + y b g 2 C 2 ijk {t) 

+ glc% k {t) + 9 2 2 ( 9 -l)ct jk (t) , 

Vb 



c} jk {t) = Re[T L ± .* fc (t)r±.(0(g* fe cot/3 + ^ fc )] , 
Cf jk (t) = M^(t)T% j (t)(Q* jk catp + S; k )] , 
C} jk (t) = Re[T L ± ; fc (t)T± y .(t)(g fcj cot/? + 5 fcj )] , 
Cf jk (t) = Re[R^ k (t)T^(t)(Q kj cotp + S kj )] , 



(A.23) 



(A.24) 
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and 



C] jk {t) = cat PMT&WT^mPjtf) + P' jk {b){ » H ta °f + Ab ))] , 

C- Jk (t) = - cotpRe[Rt ki (t)T^(t)(P jk (b) + P' jk {b){ ^^-^ ))\ , ^ 

C% k (t) = cot2/3Re[T±*(t)T±,(t)P/ fc (6)] , 
Ct jk (t)= -cot2/31te[R*(t)T^(t)P? k {b)] ■ 

The matrices and T^ 1 involve products of squark and chargino mixing matrices 

R%{t) = 0\ x XJ] x , T L %(t) = O b i2 U; 2 , T^(t) = o b tl v j2 , 

= OlVfr , T±.(&) = 0\ 2 V* 2 , T±.(fe) = 0^ 2 , (A.26) 
^) = V?i , = , T± (r) = C/ i2 . 

Finally the gluino contribution depends on Gij(t), defined by 

' H H cot/3 + 



CM = P^it) Uj^ + Pljit)^- 



m t 



g 2 cot 2 (3 — 1 , „ 



(A.27) 



The finite corrections to the /lo^fr vertex are obtained from those for the top 
vertex with the following substitutions. First interchange b t and cos (3 <-> sin /3 
everywhere in Eqs. (A. 14) - (A. 17) and in the definitions of N l , C % (and their 
L <-> i? counterparts), iV*, C* and G. Then in the substitute — Q" <-> S" 
and in the C l substitute Qj k — > 5^, Sjfc — > Qfcj- Also in Nfj k ,5 (b) and Gij(b) 
one must change the signs of the terms with tan 2 (3 — \ factors and one must 
change the overall signs of Cfj k (b). The finite corrections to the h$fT vertex are 
obtained from those for the bottom vertex by substituting b — > r, yt, mt — > 
and t — > z/ r everywhere (omitting the appropriate sums). The charged Higgs and 
gluino contributions vanish in this vertex, while the only chargino contributions 
come from wino induced contributions involving g 2 C 2 and g 2 C^. 
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We may now construct the threshold matching functions for the third genera- 
tion Yukawa couplings. Using the results of Section 2, 



\SUSY\ 

^Vt > 

, SUSY\ 



y s b M (ri = y b (fi)cos(3(l + A s y " SY ), 
y r(ri = yMcosf3(l + A s « SY ), 



(A.28) 



where 



A SUSY 



&a ~ ^(KaL + K aR + K ho ) 



(A.29) 



The explicit form for the complete matching functions is rather unwieldly however 
simplified forms can be constructed in certain limits. For example, in the low to 
intermediate tan (3 range it is appropriate to include order y\ and g| effects as the 
dominant contribution. We shall consider a particular limit in which the neutralino 
and chargino mass matrices have a definite form, determined by assuming values 
for the Higgs potential parameter /i H and the gaugino masses A4i ; 2 which are 
significantly larger than Mz- The the neutralino mass matrix has the approximate 
form 



M 2 





M x o 






v o 









o \ 



J 



(A.30) 



so that Zij = Sij, Zi3 = Z{4 = and Z;+2,j+2 = Z 2 ij, where 



Z 2 



( 



1 



1 



(A.31) 



V% \1 -1 

Due to the appearance of a negative Majorana mass eigenvalue we have = —£4 = 
sgn(/i H ) (note also = 1)- The chargino mass matrix is approximately diagonal, 





M x± 




(A.32) 



with Uij ~ Sij and V ~ diag(l, — sgn(/x H )), where sgn(/i H ) = fi H /\fi H \. Thus in this 
limit the neutralino eigenstates are the bino, neutral wino and a Dirac Higgsino 
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with masses m A4\, M. 2 and m^ m \/i H \. The chargino eigenstates 

are just the charged wino and Higgsino with masses mw ps M. 2 and m^ ± ^ \\i H \. 

In this limit, the explicit form of the vertex corrections can be obtained using 
simplified forms of the N l , N l , N l , C l , C l and C l . For the neutralino contributions 
we have 

a 



E = ^((°n) 2 *« + ^H(-(Ol 1 ) 2 + 4(Of 2 ) 2 )) , 

E = ^((°'i) 2 ^2 + |**i((Ofe) 2 + 2(Of 2 ) 2 )) , 

E 4*( t ) = ^((°ri) 2 ^2+^i(-(o[ 1 ) 2 +2(o[ 2 ) 2 )) , 

i=3,4 ^ 

Nf Jk it)= - it% lP ; k (t)(p jk (t) +p; k m llHCO l P+ ' 
Nf Jk (b) = it\ lP ; k (b) ( Pjk (b) + ^( < gta] ^ + A )) , 

*§*(r) = ~ t 2 *^ (r) fe(r) + ^(r)( ^^ + ^ )) . 



The simplified chargino contributions can be obtained using 



Cf jk (t) = -^(O b tl ) 2 6 kl 5 j2 , 



1 



y/2 

c fjk( T ) = yf 5fc1 ^ 2 , 

C} jk {t) = - sgn(^)fe cot pP' k {b) (P jk (b) + P' jk {b) + 



(A.33) 



(A.34) 



Cf jk (t)= -sgn( f i H )5 l2 cot2(3P; k (b)P;' k (b) , 
Cf jk (b) = -sgn(^)fe cot 2f3P; k (t)P;' k (t) . 

The contributions of N 2 and C* 2 are determined by multiplying the results for N 2 
and C 2 by — sgn(/i H ) cot /3 for the top case and by — sgn(/i H ) tan /3 for the bottom, 
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r cases. The other neutralino and chargino contributions vanish in this limit. The 
wavefunction renormalization contributions also simplify in an obvious way; in 
particular, the terms proportional to g 2 y a vanish in this limit. 

Below we give the dominant contributions to the matching functions and ignore 
those contributions of less than 1% to the Yukawa couplings. Some of the yi, and 
y T corrections will be suppressed in the low tan/5 region while receiving large 
enhancements for high tan/3. The converse is true for the corrections to y t . The 
gluino contributions dominate the top and bottom Yukawa corrections. For large 
tan (3 the finite vertex contributions from graphs with chargino and neutralino 
exchange are enhanced for the bottom and r cases and give large contributions. 
For high(low) tan/5 the finite parts of wavefunction renormalization contributions 
depending on yt(yb, Vr) can be important as well. The matching functions are 



167T 2 A^ sy ~ §^(-2m^(m t G2(mJ,m2,m?)^ i (t)^i(*) 
+ {n H cot (3 + A t )G 2 {m% m\, m|)(i%(t)) 2 ) 

, * (A.35) 
+ 3 Vt ( E F * , m% ± ) (O b n ) 2 + F 2 H , m% ) 

^ i i 

+ cos 2 ''{3 F 2 (m 2 H ±, m 2 )^ 

^ i ' 
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16tt 2 A^ ~ yii^-2m~ g ^ H tanP + A b )Y^G 2 (m%ml,mlXPl j (b)f 

i ' 

+Vt I ~ 2 Vh tan/? (jn t G 2 (m 2 q± , m~ , m|ji^(t).F^(t) 

+ cot /3 + A t )G 2 (m| ± , m|, m| ) (i*. (t)f 
+ \(^F 2 {ml,m\ ± ){0\ 2 ) 2 + cos 2 (5 F 2 {m 2 H± ,m 2 )) j (A.36) 

i ' 

+bl (E F2 K' rn% ± ) (Oj,) 2 + £ F 2 (m 2 bi ,m% ) 
+ sin 2 /3 F 2 (m 2 H ±,m 2 ) 



+\gl UhM 2 tan/3 ^ (c 2 (m? , m| , m^) (O^) 

^ i 

+ G 2 (m%m\ ±) m%)(0\ l ) 2 )^ , 



167r 2 A^ sy ~ \y 2 A F 2 {m% T ,m\ ± ) + ]T F 2 (m|, m| ) + sin 2 /? i^m 2 ^, 0) j 

+ ^ 2 2 ('/i^2t a n/?^(G2(m|,m| ,m|,3)(O[ 1 ) 2 (A.37) 

+ 0i «<„|U„ 

- 2t 2 .M ! (/x H tan /3 + A T ) £ G 2 (m| , m? , m? . ) (i*- (r) ) 2 ) , 



where we have used Eqs. (A. 30) - (A. 32). 
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APPENDIX B: Function Definitions 



In the appendix we define the various functions arising from parameter integrals 
in the evaluation of the finite parts of one loop diagrams. The functions Fj and G{ 
come from two and three point functions, respectively. The Fj are 

M\a + M 2 B {1 - a) 



Fi(M|,M|) = J da In 



F 2 (M A ,M B ) = 2 / daaln 



M\-M B \ A /i 2 B fx 2 
M\a + M|(l - a) 



(^^((Ml^MiM^ln^.Mlln^) ^ 
M 2 B x 



M\-M 2 B 2 ' 



F3(M|,M|) = 6 y daa(l-a)ln(^ ^ 



1 'm\{M 2 a -?,M b ) In M ^ 



(M 2 - M| )3 V"^"^ 4 "'"^ ^2 

4/ 9 9n Af|\ iM/i.l/r. 5 

-M^(Mg-3Ml).n^) + (M ^ M | )2 --. 

Both Fi and F3 are symmetric. The explicit forms of the functions for the limits 

M A > (<)M B are 

Fi(Mi,0)= la— ^-1 , 
/i 

F 2 (Mi, 0) = In ^ - i , F 2 (0, M\) = In ^ - § , (B.2) 

/i /! 

F 3 (Mi,0)= ln-^-|, 



G2 



and for M A = M B , F { = \nM\/^ 2 . The d are 



G 2 (M 2 A ,M 2 ,M 2 ) 



1 l-Q 



G 1 {M 2 A) M 2 B) M 2 ) = 2 



J da J d(3\n 





M\a + Ml/5 + M 2 (l -a-py 



-(ln^-I^ 

2 _ n /r 2\n/r2 _ T\/r2\y ni ,,2 2- 



Ml 



+ (A <-> B) + (A <-> C) - 1 , 
1 l-a 

da y dp 
o o 
1 



1 



M\a + M|/3 + Mg(l -a-P) 



M 2 



Ml 
B In B 



Mi 



M 2 - Mj, \M 2 - M 2 M 2 M 2 - M 2 M 2 



:b.3) 



We also consider the limits below 



GriMlMlM 2 ) 



2M\ / 1 M\ 

+ (Ml-M£) 
I) 



n m a 
In — f 



G 2 {Mi,M B ,M z c )= { 



M.\ Mil 

M.\ Mn » M c 

.1/,. M B > M c 
.W.i > M B ,M C 

t~ — 

(Ml M tf^ ln ^+ M^M^ M A = M B 

M A = M B ^> M c 



M\-Ml 



3 
2 



In ^4 



„ 1 „ It, A# 



(B.4) 



where we have included the particular cases that occur in the threshold corrections. 
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FIGURE CAPTIONS 



1) Vertex diagrams contributing to dominant gluino induced threshold correc- 
tions to yt and y b . 

2) Vertex diagrams involving squarks, Higgsinos and winos giving the dominant 
contributions to yj, threshold correction which are enhanced in the large tan (3 
region. 

3) Vertex diagrams involving sleptons, Higgsinos and gauginos giving the domi- 
nant contributions to y T threshold correction which are enhanced in the large 
tan (3 region. The bino exchange diagram can be important due to the large 
hypercharge of the r lepton. 

4) Variation of M GUT with gluino mass for 120 GeV < nig < 1 TeV and central 
values of sin 9w(Mz) and a(Mz)- Contours are given for the allowed range 
of a s (Mz) and two values of tan/? for fixed M t = 165.4 GeV. 

5) Variation of M GUT with M t for 130 GeV < m t (M z ) < 200 GeV and central 
values of sin 0w(Mz) and a(Mz). Contours are given for the allowed range 
of a s (Mz) and two values of tan (3 for fixed nig = 400 GeV. 

6) Variation of M Hs with V SUSY ~ .08m^ for 100 GeV < < 1 TeV and 
central values of sin 6\y(Mz) and ct(Mz). Contours are given for the allowed 
range of a s (Mz) and two values of tan (3 for fixed M t = 165.4 GeV. 

7) Variation of M„ 3 with M t for 130 GeV < m t (M z ) < 200 GeV and central 
values of sin^y^M^) and a[Mz\ Contours are given for the allowed range 
of a s (Mz) and two values of tan/3 for fixed V SUSY = 45 GeV. Note that the 
initial fall of M Hi with Mt is due to the electroweak threshold while its rise 
for large M t is due to the effect of large y t on the two loop gauge (3 functions. 

8) The dependence of the M t solution on tan/3 showing the effect of a large 
splitting of M v and M E for M b = 5.2 GeV, a s (M z ) = 0.120 and M SUSY = 
Mz- The sparticle spectrum parameters are fixed at nig = M3 = 1 TeV and 
V SUSY = 80 GeV. The typical GUT masses are then M GUT = I0 16 - 1±ai GeV 
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and M H3 = io 15 - 7±0 - 4 GeV for these solutions. We also indicate the bound 
from the nonperturbative limit on yt- 

9) The dependence of the Mt solution on tan f3 showing the effect of a splitting 
of M Hz above M GUT for M b = 4.9 GeV, a s (M z ) = 0.127 and M SUSY = M z . 
The sparticle spectrum parameters are fixed at rrig = M3 = 1 TeV and 
V SUSY = 80 GeV. The typical GUT masses are then M GUT = I0 16 - 2±ai GeV 
and M H:i = I0 16 - 8±a3 GeV for these solutions. We again indicate the bound 
from the nonperturbative limit on yt. 

10) The dependence of the Mt solution from Yukawa unification on tan/5 for the 
allowed range of a s (Mz) for M b = 4.9 GeV and M SUSY = Mz- We give 
the solutions with and without GUT scale Yukawa thresholds for the generic 
case M y 

11) The dependence of the M t solution on tan/? for bottom pole masses in the 
range 4.7 GeV < M b < 5.3 GeV for a s (M z ) = 0.120 and M SUSY = M z . 
Again, the solutions with and without GUT scale Yukawa thresholds corre- 
spond to the generic case M v = M E . The curves for M b = 4.7 and 4.8 GeV 
with GUT thresholds are cut off due to the absence of perturbative solutions. 
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